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Abstract We consider specific quantum mechanical model problems for which perturbation theory 
fails to explain physical properties like the eigenvalue spectrum even qualitatively, even if the asymp- 
totic perturbation series is augmented by resummation prescriptions to "cure" the divergence in large 
orders of perturbation theory. Generalizations of perturbation theory are necessary which include in- 
stanton configurations, characterized by nonanalytic factors exp(— a/g) where o is a constant and 
g is the coupling. In the case of one-dimensional quantum mechanical potentials with two or more 
degenerate minima, the energy levels may be represented as an infinite sum of terms each of which 
involves a certain power of a nonanalytic factor and represents itself an infinite divergent series. 
We attempt to provide a unified representation of related derivations previously found scattered in 
the Uterature. For the considered quantum mechanical problems, we discuss the derivation of the 
instanton contributions from a semi-classical calculation of the corresponding partition function in 
the path integral formalism. We also explain the relation with the corresponding WKB expansion of 
the solutions of the Schrodinger equation, or alternatively of the Fredholm determinant det{H — E) 
(and some explicit calculations that verify this correspondence). We finally recall how these conjec- 
tures naturally emerge from a leading-order summation of multi-instanton contributions to the path 
integral representation of the partition function. The same strategy could result in new conjectures 
for problems where our present understanding is more limited. 
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Chapter 1 

Introduction 



The divergence of perturbation theory in higher orders and the asymptotic character of the perturbation series 
have led to fairly extensive investigations on large-order perturbation theory [1]. Here, we are concerned with 
time-dependent, classically allowed field configurations (or classical particle trajectories that solve the Euclidean 
equations of motion) for which the contribution to the field energy (the classical Euclidean action) is finite. These 
configurations give dominant contributions the generating functionals (partition functions) in quantum field the- 
ory (quantum mechanics). The instantons manifest themselves in a particularly clear manner in the path integral 
representation from which they were originally derived. However, the analysis suffers, even in comparatively sim- 
ple cases, from a number of mathematical subtleties, and confusion may easily arise with regard to a number of 
necessary analytic continuations that assign a well-defined meaning to the resulting infinite series. Here, we are 
concerned with a number of quantum mechanical model problems in which instantons play a crucial role: in the 
cases considered, perturbation theory fails to account for the eigenvalue spectrum even qualitatively; two or more 
distinct energy levels are characterized by one and the same perturbation series. The degeneracy which persists 
after the appUcation of perturbation theory is lifted only when instanton configurations are additionally taken into 
account. The example to be considered is the quantum mechanical double-well oscillator in which the particle 
trajectory may describe an arbitrary number n of oscillations between the two minima (and we also consider gen- 
eralizations of this problem with more than two minima). The Euclidean action of two oscillations between the 
minima is found to differ slightly from the sum of two single oscillations, the difference giving rise to a modifica- 
tion of the contribution to the path integral and the partition function: this is the origin of the instanton interaction 
and the multi-instantons. The result of the investigations is that the energy eigenvalue is given as an infinite sum 
over series each of which characterizes a particular configuration in which the particle oscillates n times between 
the minima, including the modifications to the resulting Euclidean action (the n-instanton configuration). 

Let us consider a double-well oscillator with a coupling g whose inverse is proportional to the spatial separation of 
the two minima. An expansion in powers of g then corresponds to an expansion about the case of infinite separation 
of the two minima (i.e., about the case where the other minimum at x = oo can be ignored and a single minimum 
results). This cannot be assumed to be satisfactory because a situation with two degenerate minima is qualitatively 
different from a potential with a single minimum. Instanton configurations provide the necessary supplement to 
perturbation theory, which is a pure expansion in powers of g. The expansion of the eigenvalue in the instanton 
formalism is systematic and leads to a well-defined representation as a triple series in (i) the nonperturbative 
expression exp{—l/g), raised to the power n of the number of instanton oscillations, (ii) logarithms of the form 
ln{—g) which are caused by proportionality factors multiplying the instanton interaction, and (iii) powers of the 
coupling g. The logarithms ln(— g), for positive coupling g > 0, give rise to imaginary contributions which are 
compensated by imaginary parts of generalized Borel sums of lower-order instanton configurations, resulting in a 
real energy eigenvalue, as it should be. The necessary analytic continuations have given rise to some discussion 
and confusion in the past; today, we can verify the validity of this procedure by comparing the "predictions" of 
the multi-instanton expansion, including the inherent analytic continuations, with numerically determined energy 
eigenvalues of up to 180 decimal figures. 

Let us briefly expand on the connection between the generalized expansions discussed here and the theory of 
differential equations as developed by Riemann. It is weU known that the hypergeometric equation is a second-order 



5 



CHAPTER 1: INTRODUCTION 



differential equation with three singularities at ^; = 0, 1, c». The solutions of this equation, the hypergeometric 
functions, can be expressed as power series in z. If one assumes further singularities, the solutions are no longer 
expressible in terms of simple power series. Instead, they may be expressed as series of the form c (In zy 
(c constant, a € C, p integer, p > confined). The Riemann-Hilbert theorem shows that they are essentially 
characterized by the way in which they reproduce themselves under analytic continuations along closed paths in 
the complex plane. An equation with essential singularities admits solutions which are expressible as a sum of 
terms of the form c exp[g'(2) z~'^] z" (In^)^' (c constant, q a polynomial, d positive and integer, a G C,p integer, 
p > confined). Inspired by the ideas expressed in the volume [2], we call expansions of this type "resurgent 
expansion." 

Here, we discuss in detail a set of conjectures about the complete form of the perturbative expansion of the spectrum 
of a quantum Hamiltonian H, in situations where the potential has degenerate minima [2-6] . It will be shown that 
the coupling parameter g takes the formal role of the natural unit of action ft, and the systematic expansion about 
the classical problem (expansion in power of ft in the form of a WKB expansion) will be used in order to derive the 
resurgent expansions in a systematic way. Perturbative expansions are obtained by first approximating the potential 
by a harmonic potential near its minimum. They are expansions for ft ^ valid for energy eigenvalues of order ft, 
in contrast with the WKB expansion where energies are non-vanishing in this limit (WKB approximation applies 
to large quantum numbers). When the potential has degenerate minima, perturbation series can be shown to be 
non-Borel summable. Moreover, quantum tunneling generates additional contributions to eigenvalues of order 
exp(— const. / ft), which have to be added to the perturbative expansion (for a review and more detail about barrier 
penetration in the semi-classical linnit see for example [7]). 

Therefore, the determination of eigenvalues starting from their expansion for ft small is a non-trivial problem. The 
conjectures we describe here, give a systematic procedure to calculate eigenvalues, for ft finite, from expansions 
which are shown to contain powers of ft. In ft and exp(— const. / ft), i.e. resurgent expansions. Moreover, general- 
ized Bohr-Sommerfeld formulae allow us to derive the many series which appear in such formal expansions from 
only two of them, which can be extracted by suitable transformations from the corresponding WKB expansions. 
Note that the relation to the WKB expansion is not completely trivial. Indeed, the perturbative expansion corre- 
sponds (from the point of view of a semi-classical approximation) to a situation with confluent singularities and 
thus, for example, the WKB expressions for barrier penetration are not uniform when the energy goes to zero. 

These concepts will be applied to the following classes of potentials, 

• the double-well potential, 

• more general symmetric potentials with degenerate minima, 

• a potential with two equal noinima but asynometric wells, 

• a periodic-cosine potential, 

• resonances of the C(;^)-symmetric anharmonic oscillator, for negative coupling, 

• eigenvalues of the ©(z/)- symmetric anharmonic oscillator, for negative coupling but with the Hamiltonian 
endowed with nonstandard boundary conditions, 

• a special potential which has the property that the perturbative expansion of the ground-state energy vanishes 
to all orders of the coupling constant (Fokker-Planck Hamiltonians). 

Based on the extensive mathematical work carried out on the subject [2], it is perhaps rather natural to remark 
that several of these conjectures have meanwhile found a natural explanation in the framework of Ecalle's theory 
of resurgent functions [8] and have now been proven by Pham and his collaborators [9-12]. In what follows, we 
always assume that the potential is an analytic entire function. This is also the framework in which the available 
mathematical proofs exist. We first explain the conjecture in the case of the so-called double-well potential. Note 
that, in what follows, the symbol g plays the role of ft and the energy eigenvalues are measured in units of ft, a 
normalization adapted to perturbative expansions. Although several conjectures have now been proven, we believe 
that heuristic arguments based on instanton considerations are still useful because, suitably generalized, they could 
lead to new conjectures in more compUcated situations. 
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Instantons have played a rather significant role in explaining unexpected behaviour in several aspects of quantum 
field theory, for example the nontrivial structure of the quantum chromodynamic vacuum or the absence of the "su- 
perfluous" rj meson. Instantons effectively eliminate the the U (1) symmetry that would generates the rj particle, 
eUminating the problem. In the case of the one-dimensional problems studied in the current review, the instanton 

contributions can be subjected to a systematic study, including an investigation of the interaction among instantons 
which leads to a systematic expansion of the energy eigenvalues in ft. In h and cxp(— const. jK). 

This review is organized as follows. In chapter 2, we describe the conjectures. In chapter 3, we discuss the connec- 
tion with the WKB expansion as derived by different methods from Schrodinger equation. In chapter 4, we explain 
how these conjectures were suggested by a summation of the leading order multi-instanton contributions to the 
quantum partition function Tr exp(— in the path integral representation. In chapter 6, we discuss in the same 
way the periodic cosine potential, and in chapter 5 more general potentials with degenerate minima. In chapter 7, 
related problems in the example of the 0{i') quartic anharmonic oscillator are considered. Finally, the appendix 
contains additional technical details and explicit calculations verifying some implications of the conjectures at 
second WKB order. 
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Chapter 2 



Generalized Bohr-Sommerfeld 
Quantization Formulae 



2.1 Double-Well Potential 

2.1.1 Resurgent Expansion (Double- Well) 

Let us consider a quantal particle bound in a double-well potential which we write in the symmetric form 

The minima are at q = zLl/y/g. The potential is represented graphically in figure 2.1, with an indication of the 
ground state and the first excited state (g = 0.08). Further calculations profit from a translation of the coordinate 

q^q+^. (2.2) 

This translation preserves the spectrum of the Hamiltonian and leads to the quantum mechanical eigenvalue prob- 
lem 

Hil;{x) = Eijj{x) (2.3) 
in L'^{R). The Hamiltonian after the translation is 

H=-\ + V{g, q) , V{g, q)^^q^l- ^ qf . (2.4) 

For g 0, we recover the unperturbed harmonic oscillator Consequently, the unperturbed ground state in figure 
2.1 has an energy of -E = 1/2. When Rayleigh-Schrodinger perturbation theory is applied to the A^th harmonic 
oscillator state, the energy shift due to the terms of order g^ and g^ in V{g, q) [see equation (2.4)] gives rise to a 
formal power expansion in g (not ^/g), 

oo 
1=0 

For the ground state (N = 0), the first terms read 



En^J2^n19'- (2.5) 



^ 1 9 2 89 3 5013 4 

^o~2~^~2^ ~8~^ (2.6) 

The coefficients in this series grow factorially, and the series is not Borel summable (all coefficients except the 
leading term have negative sign (a complete list of the leading 300 perturbative coefficients for the ground state 
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can be found at [13]). In higher orders of perturbation theory, the numerators and denominators of the coefficients 
exhibit rather prime factors. For instance, we have 

Eo^-^^^ = -— 130116 860 668 372133 614952 623. (2.7) 

Recently, connections between number theory and physical perturbation theory have attracted considerable atten- 
tion. We will not pursue these questions any further in the current work. 

We continue the discussion by observing that for negative g, the series (2.5) and in particular (2.6) are alternating 
and indeed Borel summable. One might now expect that a suitable analytic continuation of the Borel sum from 
negative to positive g gives the desired result for the energy levels of the double- well problem. That is not the case. 
Perturbation theory fails to account for the spectrum even qualitatively. For every one eigenstate of the unperturbed 
problem (g = 0), two eigenstates emerge in the case of nonvanishing coupling g > 0. For example, the unperturbed 
ground state with an energy E = 1/2 sphts into two states with different energy for g > 0, as shown in figure 2.1. 
The new quantum number that differentiates the two states is the parity e = ±1 [see (2.10) below]. Indeed, as is 
well known, the ground state with "splits up" into two states each of which has opposite parity. The two 
levels that emerge from the ground state as the coupling g is "switched on" are separated by an energy interval of 

^ exp (-^] . (2.8) 



This term is nonperturbative in g, i.e. its formal power series in g vanishes to all orders. 
For all further discussions, we scale the coordinate q in (2.4) as 

q-^^q. (2.9) 

This scale transformation again conserves all eigenvalues. The Hamiltonian corresponding to the double-well po- 
tential can therefore be written alternatively as 

This representation shows that g takes over the formal role of h. The Hamiltonian is symmetric in the exchange 
q ^ 1 — q and thus commutes with the corresponding parity operator P, whose action on wave functions is 

Pi,{q) = ^{l-q) ^[H,P]=0. (2.11) 

The eigenfunctions of H satisfy 

H^e,N{q) = E^^N{.g)il>e,N{q) , Pi's,N{q) = £^e,N{q), (2.12) 

where e = ±1 is the parity and the quantum number N can be uniquely assigned to a given state by the requirement 
Ee,N{g) = iV + 1/2 + 0(17), i.e. by identifying the specific unperturbed harmonic oscillator eigenstate to whose 
energy eigenvalue the state energy converges in the limit g ^ 0. We have conjectured [4, 5] that the eigenvalues 
Ee,N{g) of the Hamiltonian (2.10) have a complete semi-classical resurgent expansion of the form 

^^e,.(.) = E<.^'+EG) h-iw) ^r{-^)] (2.13) 

1=0 n=i V V^9 J L \ 9J ) 

The expansion (2.13) is the sum of the perturbative expansion which is independent of the parity, 

00 

Efi9) = Y.<^3' (2.14) 
1=0 
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Figure 2.1: Double-well potential for the case g = 0.08 with the two lowest eigenval- 
ues and the unperturbed eigenstate. The energy of the unperturbed ground state is the 
energy of the (A^ = 0)-state of the harmonic oscillator ground state with E — 1/2. 
The quantum numbers of the ground state in the potential Vs{g = 0.08, q) are (+, 0) 
(positive parity and unperturbed oscillator occupation quantum number N = 0). The 
potential Vg is defined in (2.1). The energy eigenvalue corresponding to the state 
(+,0) is -E+.o = 0.317 851364 6.... The first excited state (—,0) has an energy 
of £'- .0 = 0.566 114 759 4 . . . The arithmetic mean of the energies of the two states 
(±, 0) is 0.441 983 062 . . . and thus different from the value 1/2. This is a conse- 
quence of the two-instanton (more general: even-instanton) shift of the energies. 



and of the n-instanton contributions, each of which is given by 

Here, (g) is the perturbation series whose coefficients e'"j^\ only depend on the quantum number N (and the 
order of perturbation theory I), but not on the parity e [see also (2.5)]. The index n gives the order of the "instanton 
expansion". One-instanton effects are characterized by a single power of the nonperturbative factor exp(— l/6.g), 
two-instanton effects are multiplied by cxp(— l/3g), etc. Actually, as shown below in chapters 4.3 and 4.4, the one- 
instanton effect E^^^^{g) involves a summation over n-instanton configurations, however neglecting the instanton 
interaction. The two-instanton shift E]^ ^{g) and higher-order corrections involve the instanton interaction. 

The logarithms of the form h\{—2/g) also enter into the formalism, their maximum power (for given n) is n — 1. 
Thus, the two-instanton effect carries one and only one logarithm etc. Each of the infinite series 

oo 

'^eN.nkig'' (2.16) 

1=0 

is characterized by the state quantum number N, the instanton order n and the power of the logarithm k. The series 
X) s N,nkig^ in powers of g are not Borel summable for .g > and have to be summed for g negative. The values for 
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g positive are then obtained by analytic continuation, consistently with the determination of \n{—g). In the analytic 
continuation from g negative to g positive, the Borel sums become complex with imaginary parts exponentially 
smaller by about a factor exp(— 1 /3g) than the real parts. These imaginary contributions are canceled by the imag- 
inary parts coming from the function ln(— 2/g'). For instance, the imaginary part incurred by analytic continuation 
of the perturbation series [equation (2.14) below] is canceled by the "explicit" imaginary part originating from the 
two-instanton effect [the term with n = 2 in (2.13)]. 

The instanton terms E^^^{g) are parity-dependent for odd n. In [14], a different notation was used for the coeffi- 
cients that enter into the expansion (2.13). When comparing formulas, observe that in [14], e denotes the coefficient, 
whereas here, the sUghtiy different symbol £ = ± denotes the parity, and we have 

The e-coefficients allow for a more compact notation. Their introduction is based on the observation that the odd- 
instanton coefficients have different sign for states with different parity, whereas the even-instanton contributions 
have the same sign for states with opposite parity. Thus, odd instantons contribute to the separation of energy levels 
with the same N, but opposite parity, whereas the energy shift due to even instantons is the same for both levels 
with the same N, irrespective of the parity. To give an example, we quote here from [14], without proof, several 
known coefficients for the one- and two-instanton contributions, from which the parity/sign pattern is apparent: 

,(0,+) _ ^(0 -) _ , ^(0,+) _ (0 - 
^100 ~ ^100 ~ ^ ' ^101 ~ ^101 



^1 



6299 



-101 ~ ''loi ~ 288 ' ~ ~ 

jo,+) _ (o,-)__53 (0,+) _ (0,-) _ _ 

^211 ^ ^211 — g > ^212 — ^212 — 

J0,+) _ (0 -) _ (0,+) _ (0,-) _ 23 

^200 ~ ^200 ~ 7 I f 201 — ^201 ~ 



2 



^(0,+) _ ^(0,-) _ 13 _ 1277 

^202 — ^202 — j2 ^ ' 



12 ' 


(2.18a) 


1, 


(2.18b) 


1277 

72 ' 


(2.18c) 


53 


(2.18d) 
(2.18e) 



Here, 7 = 0.57221 ... is Euler's constant. Utilizing the e-coefficients which enter into equation (2.13), these 
results can be summarized in a much more compact form as 

71 6299 
Co, 100 = 1, 60,101 = -Y^, 60,101 = — (2.19a) 

53 1277 

60,210 = 1, 60,211 = ^, 60,212 = =7^, (2.19b) 

6 72 
23 53 13 1277 
60,200 = 1, 60,201 = — ^ -^1, 60,202 = T^^' (2.19c) 

The coefficient eo,ioo gives the one-instanton separation of the levels with quantum numbers (±, 0) which can be 
approximated as 

£:_,o - E+fl « ^ exp ( - ^ ) eo,ioo, (2.20) 

^/^^ \ QgJ 

in agreement with (2.8). The well-known, classic treatment of the double-well problem is based on the "matched 
addition" or "matched subtraction" of the two solutions that characterize the two solutions in each single well. Only 
the leading term in the quantity exp(— 1/65) can be recovered in this way [15]. The advantage of the expansion 
(2.13) is that it generalizes the classic treatment to all orders in exp(— 1/65), thereby giving a clear interpretation 
of the way in which the many-instanton contributions manifest themselves in a generalized resurgent expansion 
in higher orders in exp(— 1/65). Also, the expansion (2.13) gives a systematic understanding of the error incurred 
upon terminating the instanton expansion at a given order. 

2.1.2 Quantization Condition (Double- Well) 

We have conjectured [16] that all series that enter into the generaUzed resurgent expansion (2.13) can be obtained 
by a small-g expansion of a generalized Bohr-Sommerfeld quantization formula. In the case of the double-well 
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potential, this quantization condition reads 



27r 



1 



exp 



: £1, 



or alternatively 



Ba^{E,g) 



1 -1 



+ 



£1 



\/27f V 9/ 



exp 



AdAE,9) 



0, 



(2.21) 



(2.22) 



(£ = ±1) where the A and the -B-functions enjoy symmetry properties and can be expressed as formal series in 
their arguments, 



BdAE,9) 
Ady,{E,g) 



-B^^i-E, -g) = E + J29'' bk+i{E), 



k=l 



^ oo 

-Aa^{-E, -9) = ^ + ^9^ ak+i{E) . 



(2.23a) 



(2.23b) 



fe=i 



The coefficients bk (E) and ak (E) are odd or even polynomials in of degree k, as necessitated by the symmetry 
properties. The three first orders, for example, are 



BdAE,9) = E + g (^3E^ + ^^+g^ (^35E^ + ^E^+0{g^), 
AdAE,9) = ^^"'+5 + 

V (227 E' + ^E^+0{g'). 



(2.24a) 



(2.24b) 



The function B^^jiE, g) has been inferred from the perturbative expansion, and indeed, the perturbative quantiza- 
tion condition reads 

Bi^{E,g)=N+\. (2.25) 

Here, is a noimegative integer. The perturbative expansion for an energy level is obtained by solving this equa- 
tion for E. The function Adv,{E, g) had initially been determined at this order by a combination of analytic and 
numerical calculations. 

The formula (2.13) is an exact expansion for the energy eigenvalues of the double-well anharmonic oscillator at 
finite g, which can be inferred by a systematic expansion of the quantization condition (2.22) for small g. The 
generaUzed resurgent expansion (2.13) can be derived either by a systematic expansion of the path integral repre- 
senting the partition function about nontrivial saddle points of the Euclidean action, or by a systematic investigation 
of the WKB expansion. Both of these methods will be discussed in the sequel. 

We illustrate the statements made by considering some example calculations. Specifically, we start by investigating 
approximations to the energies of the states (+, 0) and (— , 0) which emerge from the unperturbed oscillator ground 
state. The resurgent expansion is explicitly given in (2.13), with the first perturbative terms provided in (2.6) and 
the leading instanton coefficients (2.19). Neglecting higher-order perturbative (in g) and higher-order instanton 
effects [in exp (—1/65)], the energies of the two lowest-lying states are given by 



E, 



'e,0 



exp 



1 

65. 



Here, £ = ± is again the parity. We define the abbreviations 

1 



exp 



1 



and 



X{9) = In 



(2.26) 



(2.27) 



(2.28) 
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An expansion of the quantization conditions (2.21) and (2.22) for small ^{g) [note that a small positive g implies 
a small ^{g)] then involves an evaluation of the expression 



Sdw(£^e,o,5) «r(ee(r7))«e 



1 





1 


exp 


69 



(2.29) 



The quantization condition (2.21) then takes the form 



1 




fWg 




\exp 


1 

6s 



1/2 



exp 



1 1 

2 3^ 



ei. 



(2.30) 



This equation is approximately fulfilled provided that the analytic continuation of the square root is carried out 
such that for (/ > 



1/2 



1/2 



(2.31) 



A systematic expansion of the quantization conditions (2.21) or ahernatively (2.22) in powers of x(.9)^ '^(ff) and 
g then leads in a natural way to the generalized resurgent expansion (2.13). We recall that some of the known 
coefficients for the ground state have been given in the form of the t~^^:{^^ (2.18) or in the form of the bn nki- 
coefficients (2. 19). We now turn to a discussion of states with general N. The treatment is simplified if we restrict 
the discussion to the dominant effects in each order of the instanton expansion, i.e. if we restrict the discussion to 
the terms with Z = in the resurgent expansion (2.13). We employ the approximation 



in (2.23a) and obtain the following approximative quantization condition 



-l/6g 



£1 



COS ttE 



-l/6g 



: ei ■ 



2lT 



(2.32) 



(2.33) 



We show in chapter 4 that, from the point of view of the path integral representation, the successive terms corre- 
spond to multi- instanton contributions, i.e. to contributions which are generated by the modification of the instan- 
ton action due to the overlap region which results when "gluing together" two (or more) instanton solutions to a 
multi-instanton configuration. 



For example, the term n = 1 in (2.13), which is also the one-instanton contribution, is 



e''^\9) = 



e 



JV+1/2 ^_i/63 



(1 + 0(5)). 



This result generalizes the coefficient eo,ioo = 1 (2.19) to the case of arbitrary N, 

1 



ejv,ioo = 



NV 



The term n = 2, which is the two-instanton contribution, is 



E'S\9) = 



[N\Y 



9/ 



2A^+i e-i/3s 



27r 



[\n{-2/g)-^{N + l) + 0{g\ng)], 



(2.34) 



(2.35) 



(2.36) 



where if) is the logarithmic derivative of the F-function. This result generalizes the coefficient eo,2io = 1 (2.19) to 
the case of arbitrary N, 

1 ^ 



ejv,2io 



and the result 69,200 = 7 to 



ejv,200 



N\ 



(2.37) 



(2.38) 



13 



CHAPTER 2: GENERALIZED BOHR-SOMMERFELD QUANTIZATION FORMULAE 

2.1 Double-Well Potential 



More generally, it can be easily be verified upon inspection of (2.13) that the n-instanton contribution has at leading 
order the form 

E^;;\9)=[^^j (-^^j [P:{H-9m)+0{9i\ngr')], (2.39) 

in which {a) is a polynomial of degree n—1 (not the Legendre polynomial). For example, for A?^ = one finds 

P°(a) = l, P°(a)=a + 7, i^3° W = ^ + if + ^ , (2.40) 
in which 7 is Euler's constant: 7 = — V'(l) = 0.577215.... 

2.1.3 Large-Order Behaviour (Double- Well) 

The perturbation series fails to describe the energy levels at the one-instanton level. However, the mean of the 
energies of the levels {+,N) and (— , iV) may be described to greater accuracy by the perturbation series than 
each of the two opposite-parity energy levels, because the one-instanton effect has opposite sign for opposite- 
parity state. The perturbation series, which is nonalternating, may be resummed using standard techniques; the 
only difficulty which persists is the singularity along the integration axis in the Laplace-Borel integral. However, 
one may easily define a generalized Borel transform. The calculation of such a generalized transform is rather 
trivial and discussed for example in [17]. In particular, the Borel sum naturally develops an imaginary part for 
g > which represents a fundamental ambiguity. In the case of the quantum electrodynamic effective action, the 
imaginary part which describes pair production emerges from an appropriate analytic continuation. In the case 
of the double-well potential, the energy values are real, and there are compensating imaginary contributions from 
higher-order instanton effects which compensate the imaginary part incurred by the Borel sum of the nonalternating 
perturbation series. 

The mean of the energies of the levels N) and (— , N) is described by the perturbation series up to a connmon 

shift received by both levels (±, N) via the fwo-instanton effect. Thus, it is up to the fvvo-instanton level that 
the perturbation series has a physical meaning: indeed, perturbation theory describes approximately the mean of 
the energies of the levels with the same N, but opposite parity. The mean of the energy of these levels, to better 
accuracy, is then given as the sum of the real part of the (complex) Borel sum of the perturbation series, and of 
the two-instanton effect. Thus, at the level of the two-instanton effect, we expect the asymptotic perturbation series 
to go through its minimal term which, for a factorially divergent series, is of the same order-of-magnitude as the 
imaginary part of its Borel sum and also of the same order-of-magnitude as the difference of the real part of the 
generalized Borel sum of the perturbation series and the result obtained by optimal truncation of that same series. 

Thus, after an analytic continuation of the series 'YlaLa^N^nki 9^ in (2.13) from g negative to g positive, two 
things happen: the Borel sums become complex with an imaginary part exponentially smaller by about a factor 
exp(— 1/33) than the real part. Simultaneously, the function h\{—2/g) that enters into the two-instanton contribu- 
tion in equation (2.13) also becomes complex and gets an imaginary part ±i7r. Since the sum of all contributions 
is real, the imaginary parts must cancel. This property relates, for example, the non-perturbative imaginary part of 
the Borel sum of the perturbation series to the perturbative imaginary part of the two-instanton contribution. We 
have identified (in the sense of a formal power series) the quantity (2.14) 

00 

as the perturbative expansion. We now speciahze to the {N = 0)-states with positive and negative parity. The 
relation to the two-instanton contribution as given by the polynomial P2 implicitiy defined in (2.39) is 

Im4°)(ff)^^^-^e-i/3«Im [P^{\n{-g 12))] = -h-'"' ■ (2-42) 

An introductory discussion of dispersion relations and imaginary parts generated by the summation of nonalternat- 
ing factorially divergent series is given in appendix D, especially in chapter D.l. The coefficients E^^^^. are related 



14 
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to the imaginary part by a Cauchy integral [see equation (D.5)]: 



£(0) - 1 r 



Im 



(2.43) 



For A; — > oo, the integral is dominated by small g values. From the asymptotic estimate (2.42) of Im E^'^^ for 
5^0, one then derives the large order behaviour of the perturbative expansion [1]: 



1 



(2.44) 



Let us us now consider the perturbative expansion about the one-instanton shift in the sense of equation (2.15), 



with the convention -E^^^ = ejv,iofe 



(2.45) 



fe=o 



One may derive, from the imaginary part of , the large-order behaviour of this expansion, by expressing that 
the imaginary part of e'q'' (g) and E^^ (g) must cancel at leading order: 



Im4i)(ff)--£ 



-l/6g 



Im [P3(ln(-ff/2))] . 



(2.46) 



The coefficients E'f}^ again are given by a dispersion integral: 
Using then equations (2.40) and (2.46), one finds 



(2.47) 



If 



3 ( In - + 7 ) e 



9 



k+2 ■ 



(2.48) 



At leading order for k large, g can be replaced by its saddle point value l/S/c in In g and, finally, one obtains 



(I) 3'=+^ 
^kfi = eo.iofe = k\ 



ln(6fc) + 7 + 



(2.49) 



Both results (2.44) and (2.49) have been checked against the numerical behaviour of the corresponding series for 
which 300 terms can easily be calculated [13, 14]. 



2.1.4 The Function A (Double-Well) 



A specific function A of the energy levels of the double-well potential, expressed as a function of the couphng 
g, has been the subject of rather intensive investigations in the past [5]. Here, we discuss the motivation for the 
definition of A, as well as its evaluation and the comparison to the instanton expansion. 

We define the quantity B ^E^^ (g)| as the real part of the Borel sum of the the perturbation series (2.14), 



B 



{<(.)} = 



Re 



Borel sum of E^^'^ {g) = ^ e'^] g' 



1=0 



(2.50) 



In [17], some explicit examples are given for the calculation of Borel sums in cases where the perturbation series 
is nonalternating (as in our case). In the same sense, we denote by B (g)| the real part of the Borel sum of 
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the the perturbation series about one instanton in the notation (2.15), evaluated for the state with negative parity 
(so that the result, for > 0, is positive): 



B 



Re 



Borelsumof ^ eN,wig^ 



1=0 



(2.51) 



An illustrative example is provided by equation (8.13) below, where the case N = is considered and the first 
terms {I = 0, . . . , 8) are given explicitly. The quantity B ^E^\g)^ is of importance for the entire discussion in 
chapter 8.9 below. 



We now consider the ratio 



A(5) 



E-n)-B 



(2.52) 



(£;+,o-S-,o) [ln(2.9-i)+7] 

In this expression -E+,o and -E_,o are the two lowest eigenvalues of H. In the sum (-E+,o + -£'-,o) the contribu- 
tions corresponding to an odd number of instantons cancel. Therefore, the numerator is dominated for g small 
by the real part of the two-instanton contribution proportional to Re P2 . Note that, if the perturbation series were 
able to describe the mean of the energies -E+,o and -E_,o exactly, then the numerator would vanish. The difference 
— E- ) , as we know, is dominated by the one-instanton contribution which in turn is characterized by a nonper- 
turbative factor exp[— l/6.g]. Squaring the energy difference, we obtain an expression of the order of exp[— l/Sf/] 
which is characteristic of the two-instanton effect. Forming the ratio of the two-instanton effect, divided by the 
squared one-instanton, we obtain an expression of order one (the additional logarithm appears because of the 
nonvanishing eg. 210 coefficient). Using the expansion (2.13) together with the expUcit results for the coefficients 
(2.19), one verifies easily that A{g) indeed tends to unity as g ^ 0. 

Moreover, performing an expansion in powers of g and inverse powers of\n{2/g) and keeping only the first few 
terms in {l/ln{2/g)} in each term in the g-expansion, one finds [14] 



A(g) ~ 1 + 33- 



23 



135 



2 ln(2/g) 
9' 



7 



+ ■ 



7 



- O 



ln(2/s) 



ln(2/s) ln2(2/s) Vln--'(2/g 
7 , 7 



ln(2/g) ln2(2/9) 



+ 



(2.53) 



The higher-order corrections, which are only logarithmically suppressed with respect to the leading terms l + 3g, 
change the numerical values quite significantly, even for small g. 



Table 2.1: The ratio A{g) as a function of g. 



coupling g 


0.005 


0.006 


0.007 


0.008 


0.009 


0.010 


A(g) num. 


1.0063(5) 


1.0075(5) 


1.00832(5) 


1.00919(5) 


1.00998(5) 


1.01078(5) 


A(g) asymp. 


1.00640 


1.00739 


1.00832 


1.00919 


1.01001 


1.01078 



Table 2.1 displays numerical results, obtained by solving the Schrodinger equation in a range of small values of g 
for which the evaluation is still reasonably precise; there are in agreement with the first few asymptotic terms up 
to numerical precision [14]. 

Of course, for larger values, significant deviations from the leading asymptotics must be expected due to higher- 
order effects; these are indeed observed as seen in figure 2.2. For example, at g = 0.1 the numerically determined 
value reads A(O.l) = 0.87684(1) whereas the first asymptotic terms sum up to a numerical value of 0.86029. 
Further calculations are described in chapter 8, where also reference values for A{g) are given. 
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Figure 2.2: Double-well potential: Comparison of numerical data obtained for the 
function A{g) with the sum of the terms up to the order of of its asymptotic ex- 
pansion for g small [see equation (2.53)], solid line. If we express both the numerator 
and the denominator of (2.52) as a power series in g [with the help of the instanton 
coefficients (2.18)] and keep all terms up to the order of g^, the dashed curve results. 
This latter approach is also used for the curves in figures 8.1,8.2, and 8.3 below. 



2.2 Other Potentials 
2.2.1 Symmetric Wells 

We consider a slight generalization of (2.21), i.e. a potential with degenerate minima atq — and q — qo, and 

V{q) = ^u;'q^+0{q'), V {q) = ^ cj^q - qof + O {{q - qo)') . (2.54) 

The potential V is supposed to fulfill the symmetry condition V{q) = V{qo — q). Without loss of generality, one 
may set 

Lu = 1 (2.55) 
in (2.54), by a simple scaling argument. Indeed, a Hamiltonian [see equation (2.10a)] 



H{g,u^) 



satisfies 



H{g,u:) = uH(i-,l) , 



(2.56) 



(2.57) 



which maps the spectrum of the general case with arbitrary lo onto a seemingly more special case, with cj = 1 and 
a redefined coupling g. The constant C is given as [see also equation (E.lOb) below]. 



C = 9o exp 



1 



1 



1 



V2V%) 1 <io-q 



(2.58) 
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As in the case of the double-well potential, the function Asyja{E, g), which described the instanton-related effects, 
has no term of order [cf. equation (2.23b)] 



{E,g)=a/g + 0{g). 



The constant a is given by [see equation (E.lOa) below] 



2 rdq^/m^. 
Jo 



(2.59) 



(2.60) 



There is, just as in the case of the double-well problem, degeneracy on the level of perturbation theory which is 
Ufted by instantons. The quantization condition is conjectured to be 



1 



27r V2 



2(7 

B,y^{E,g) ) (-y 



-Bsym(£',5) 



exp 



^sym 



{E,9) 



■ £1, 



(2.61) 



with e denoting the parity. The function BsymiE, g) ^ E + 0{g) is related to the perturbative expansion about 
each of the two symmetric wells, as given in (3.3 1) below and explained in detail in chapter 3.2.5. The perturbative 
quantization condition reads 

B,y^{E,g)=N+l (2.62) 

with integer positive N. Explicit calculations of the first few terms in the expansion of BsymiE, g) in g are carried 
out in appendix B. 1, and relevant formulas can be found in equation (B. 11). For Agym {E, g), see appendix F.5. 

It is important to remark that the A and B-functions (2.89a) and (2.89b) are determined by the formal properties 
of the potential alone. Specifically, these functions may be inferred, for a given potential, via an evaluation of the 
contour integrals of the WKB expansion [left-hand side of equation (3.57)]. We distinguish here several different 
A and B functions by assigning indices to these symbols, depending on the potential in question. However, it is 
necessary to remember that the potential alone determine the instanton-related j4-function and the perturbative 
i?-function are determined by the potential alone [which we always use in the convention (2.10a)]. 



2.2.2 Asymmetric Wells 

In the case of a potential with two asynmietric wells, we have to modify the quantization condition (2.21). Specif- 
ically, we will discuss the case of a potential with two asymmetric wells about which one may expand 

V{q) = \u;lq^ +0{q^), V {q) = \ ujl(q ^ q^f + O {{q - q^f) . (2.63) 

Of course, the asyrmnetry manifests itself in the relation V{q) ^ V{qQ — q). The generalized Bohr-Sommerfeld 
quantization formula takes the form [see also equations (3.57), (5.1), (5.3), (5.8), (5.16), (7.20b), and (F.35) below] 



T{\-B,{E,g))V{\-B^{E,g)) 

1 / 2Ci\^''-^'^^ / 2C2\'^^^^'®^ 



2-K 



/_2^\ /_2^\ ^-A(,,E)^Q^ (2.64) 



where Bi {E, g) and B2{E, g) are determined by the perturbative expansions (see chapter 3.2.5) around each of 
the two minima of the potential 

B^{E,g) = E/uj^+0{g), B2{E,g) = E/uj^ + 0{g) , (2.65) 

The perturbative quantization conditions are 

B,{E,g) = Ni + ^, B2{E,g)=N2 + \, (2.66) 

with nonnegative integer Ni and A^2- In equation (2.64), the quantities Ci and C2 are numerical constants, adjusted 
in such a way that A{E, g) has no term of order g°: 

A{E,g)=a/g + 0{g). (2.67) 
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Again, a is given by 

rQo 

= 2 / dq ^/2V{q) . (2.68) 
Jo 

Equation (2.64) has two series of energy eigenvalues, 

En={N+^) LJi+Oig), En = {N+^) c02 + O{g) (2.69) 

close for g ^ to the poles of the two F-functions. The same expression contains the instanton contributions to 

the two different sets of eigenvalues. 

The quantization condition (2.64) is qualitatively different from (2.21) and (2.61) because parity is explicity broken 
in the asynmietric case. One verifies that multi-instanton contributions are singular for a; = 1. But, if one takes the 
symmetric limit, one obtains a form consistent with the product of the two equations (2.61). 

Without loss of generality, one may set 

Wi = 1 , ij02= ijJ (2.70) 
in (2.63), by a simple scaling argument. Indeed, a Hamiltonian [see equation (2.10a)] 

H{g, a;2) = - 1 ( + i Q ^2 ^2 ^ 1 ^2 _ ^^^2^ (2.71) 

satisfies 

/f(5,a;i,W2)=wi-fff— ,1,— ) =H{g,l,uj), (2.72) 

a relation that maps the spectrum of the general case with arbitrary g, wi, and W2 onto the more special case g, 
uji = 1,LJ2= ij^. One may than redefine g ^ g and a; — > w. The constant C, defined as [see equation (5.3) below]. 



2w 

l+UJ 



n° ( I 1 1 

^0 



(2.73) 



/o \y/W(^ q w{qo-q)^ 
tends to (2.58) in the limit u) ^ 1. The constants Ci and C2 are then given by [see equation (5.8) below]: 

Ci=C^, C2 = C^/uj. {2.1 A) 

These enter into the quantization condition (2.64). We note also that the quantization condition (2.64) may be 
derived, to leading order in gi, by an exact evaluation of the multi-instanton contribution to the path integral [see 
(5.16)]. 

We now consider the imaginary part of the leading-instanton contribution and large-order estimates. The expression 
(2.64) can be used to determine the large order behaviour of perturbation theory by calculating the imaginary part 
of the leading instanton contribution and writing a dispersion integral. Setting wi = 1, a;2 = w, for the energy 
^Nig) = N + ^ + 0{g) one finds the imaginary part 

ImEj^ig) - if^^-(JV+i/2)(i+i/c.)g-a/s (2.75) 

with [equation (2.73)] 

= izlL^^-(iv+i/2)/.. (2^)^^+'/'^^'+'/"' (2.76) 
2iT Ni 

X sin [jt{N + i)(l + 1/uj)] F [i - (Af + i)/a;] . (2.77) 

From lm.EN{g), one infers that the coefficients E^j^\ of the perturbative expansion of E^ig) behave, for order 
fc — > 00, Uke 

p{0) ^ r[fc + (iV + l/2)(l + l/a;)] , 

Note that this expression, in contrast to the instanton contribution to the real part, is uniform in the limit w = 1 in 
which the result (2.44) is recovered. 
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2.2.3 Periodic Cosine Potential 

We first consider the periodic cosine potential. This example differs from the preceding ones because the potential 
is still an entire function but no longer a polynomial. On the other hand the periodicity of the potential simplifies 
the analysis, because it allows classifying eigenfunctions according to their behaviour under a translation of one 
period: 

V^((7 + T)=e'^V^(g), (2.79) 

where T is the period. For the cosine potential ^(1 — 00849) (and thus T = 7r/2), the conjecture then takes the 
form 

2x-Bpc(E,s) eAp.(B,g)/2 / 2nBpc(E,s) ^-A^^(E,g)/2 _ 3 COS 

r[i-Spe(£,ff)]^Vy; T[^ + Bp,{E,g)] ~ ' ^^'^^^ 

The first few terms of the perturbative expansions of the functions and -Bpc for the periodic potential are 

Bp,{E,g) = E + g (^E^ + ^'j+g^ (^3E^ + ^E^+0{g^), (2.81a) 

Ap,{E,g) = g~^+g(^3E^ + ^'j +g' (^UE^ + ^E^ +0{g'') , (2.81b) 

where again A has been initially obtained by a combination of analytic and numerical techniques. The perturbative 
quantization condition is 

Bp4E,g)=N+^. (2.82) 

The conjecture has been verified numerically up to four-instanton order with decreasing precision. The theory of 
resurgence has allowed to also investigate this potential, as well as general trigonometric potentials [12]. 



2.2.4 (!?(z^)-Anharmonic Oscillator and Fokker-Planck Equation 

We consider now the 0{u)-symmetric anharmonic oscillator with the Hamiltonian 

/f = -iv2 + iq2+g(q2)^ (2.83) 

where q belongs to W . The analytic continuations of the eigenvalues to <? < are complex, the imaginary part 
determining the decay rate by barrier penetration of states concentrated initially in the well at q = 0. Moreover, 
the behaviour of the imaginary part for 51 — > 0_ is directiy related to the large-order behaviour of the perturbative 
expansion. 

The eigenfunctions of H can be classified according to the irreducible representations of the orthogonal 0{v) 
group. We work in IR'^; at fixed angular momentum I, the Hamiltonian (2.83), expressed in terms of the radial 
variable r = |q|, takes the form 

i/^ay iv-id^ii(i + v-2) i^ 4 



After the change of variables 
one finds 

with 



Hi^Hi= rt^-i^/^ Hi r^^-^^l^ , (2.85) 
B.--\{^WV^^\r^^9r^ (2.86) 



2 Vdr 



j = Z + Z//2 - 1 . (2.87) 



The Hamiltonian Hi is one-dimensional with a non-polynomial potential. Eigenvalues depend on the parameter 
j = (but with the constraint > u/2— 1). Note also the formal symmetry j ^ —j. 
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Writing the perturbative quantization condition as 

B^{E,g,j) = l+j + 2N, N>0, (2.88) 
one finds [see also equations (2.24a), (7.16), and (7.48)] 

BAE,g,j)=E-^ [3E^-f + l]g+^ [7 E^ + [5 - 3 f) E] + O {g') . (2.89a) 
The first few terms of the expansion of the function A^, are [see also (7.49)]: 

A4E,g,j) = --g-^+l^-j^----E^j g 

+1 (227 E^ - 77 f + 187) E g^ + O (g^) . (2.89b) 
8 

The A and B functions are determined exclusively by the semi-classical expansion of the solutions of the spectral 
equation of the Hamiltonian (2.84). 

In order to write down a quantization condition, it is necessary to specify boundary conditions for the Hamiltonian. 
In the current chapter, we will take advantage of a well-known connection between the O (i/)-anharmonic oscillator 
at negative (!) coupUng on the one hand to the double- well potential at positive coupling with a symmetry-breaking 
term on the other hand [18, 19]. (This result, initially conjectured on the basis of numerical evidence, was first 
proven later by path integral manipulations [20] or recursion formulae [21], and later generahzed to arbitrary 
3 [19].) 

We briefly recaU here the central statement regarding the connection between "broken-double-well and the 0{p)- 
symmetric quartic potential. Indeed, as shown below in chapter 7.2.1, the spectral equation implied by the Hamil- 
tonian 

I f dV I 1^-1 d 1 1(1 + 1^-2) I2 4 



may be reformulated under appropriate substitutions into 

1 / 1 \ ^ E 

The spectrum of this equation is bound from below and leads naturally to a self-adjoint extension of the 0(1/) 
anharmonic oscillator at negative coupling. This spectrum is not the same as the resonances implied for negative 
coupling if the Hamiltonian is endowed with boundary conditions that imply a vanishing of the wave function 
in complex directions of the parameters, which would lead to resonances discussed below in chapter 2.2.5. As 
is evident from (2.90a) and (2.90b), the Hamiltonian Hi for g < 0, with proper boundary conditions, thus has a 
self-adjoint extension, which in turn is (almost) equivalent to the double- well potential, but with a linear symmetry- 
breaking term. 

In the case j = the linear symmetry-breaking term in (2.90b) vanishes. Setting q = p + 1/2, equation (2.90b) 
becomes identical to the spectral equation of the double-well potential. This consideration leads immediately to 
the foUowing correspondence 

B,{E,-g,j = 0) = 2BdUE/2,g), (2.91a) 
A4E,-g,j = 0) = Aa^{E/2,g). (2.91b) 

The correspondence between the j = 0-component of the 0(i/)-anharmonic oscillator and the (symmetric) double- 
well potential [18, 19], thus finds a natural extension to instanton effects described by the function Ai,{E, —g,j = 
0) = Adv,{E/2, g). The perturbation series implied by the condition Bd^{E/2, g) = N + \ leads to a non-Borel 
summable series for £ as a function of g; the imaginary parts are compensated by instantons. 
In view of (2.89a), (2.89b) and (2.90b), we conjecture for the self-adjoint extension of the Hamiltonian (2.84) at 
negative coupling —g (with 5 > 0) the quantization condition 

1 



r [i(l + j - B,{E, -g,j))\ r [i(l - j - B,{E, -g.j))] 



gj 2-K 
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In leading order in g, this equation may be derived via instanton calculus [see equation (7.46)]: 

1 / 2\''e-i/3f 



T{^il+j-E))r{lil-j-E)) V 9 J 277 



. (2.93) 



We now turn our attention to an important special case. Setting j = — 1 in (2.90b) and shifting E 2E, we obtain 
the Fokker-Planck Hamiltonian 



with the associated eigenvalue problem 

HFPipip) = Eip{p), (2.94b) 

which is known to have a peculiar property: the perturbative expansion of the ground-state energy vanishes to all 
orders in the coupling. Of course, the symmetry j — j, which is present in (2.86), imphes that the the spectrum 
of (2.94) is invariant under a change of the sign of the linear symmetry-breaking term. 

{Remark. ) With the shift p = q — l/2m (2.94), we may rewrite the Fokker-Planck Hamiltonian as 



H - ^ ^ + ^ 



{1-qf + gq 



\ . (2.95) 



We may thus identify the Fokker-Planck Hamiltonian as originating from the double-well Hamiltonian (2.10b), 
with an additional linear symmetry-breaking term of relative order g, and an additional global shift of — 1 /2 which 

compensates the ground-state energy of the harmonic oscillator 

In view of the correspondence implied by (2.90) and the necessary shift E 2E in going from (2.90b) to 
(2.94), we conjecture the following secular equation for the Fokker-Planck Hamiltonian as a special case of the 
anharmonic oscillator (2.92), 

/_2y--(-^^)exp(-^.p(^,,))^^_ (2.96) 



T{-BFp{E,g))T{l-BFp{E,g)) \ gj 2-k 
The perturbative quantization condition is 

BFp{E,g)=N (2.97) 
where A?" is a nonnegative integer. The functions Bpp {E,g) and ^fp {E,g) aie given by 

B4E,-g,j = ±l) = 2Bpp(^^,g^, (2.98a) 

A^{E,-g,j = ±l) = AFp(^|,,g), (2.98b) 

where either sign of j leads to the same result. The first terms read 

Bpp {E,g) = E + 3E^g+ (35 E^ + ^E^ g^ + 0{g^) , (2.99a) 

AFp{E,g) = 1-+(^17E^ + ^^ g+ (^227 E^' + ^E^ g^ + 0{g^). (2.99b) 

Indeed, the quantization condition (2.96) may be derived, at leading order in g, using instanton calculus [see 
equation (7.29) witii C = 1]: 
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2.2.5 O(zy)-Anharmonic Oscillator: Resonances 

In the current chapter, we investigate the quantization condition for the resonances of the C'(^)-anharmonic oscil- 
lator at negative coupling. These may be found by analytic continuation of the eigenvalues to negative coupling. 
The notation is the same as in chapter 2.2.4. 

First, we briefly investigate the special case j = at g > 0. The correspondence (2.91) implies the perturbative 
condition B^{E,g,j = 0) = 2Bdw(E/2, —g) = 2N + 1 for eigenvalues of the e'(i^)-anharmonic oscillator at 
positive coupling. 

The generalization to arbitrary j at g > implies a Borel-summable, alternating series for E as a function of g. 
Indeed, the perturbative quantization condition (2.88) becomes exact in this case, 

B4E,g,j) = l+j+2N, N>0. (2.101) 



We now turn to the resonances for arbitrary j atg < 0. We recall that the A^, and i?^-functions (2.89a) and (2.89b) 
are determined by the formal properties of the potential alone. However, the actual quantization condition may also 
depend on the boundary conditions. The secular equation for resonances at (; < is conjectured to be 



le 



-A.{E,g) 



exp 

T[\ {j + l-B,{E,g))\ 
■V[\ {j + l + B,{E,g))\ 



ITT 



j + l+B,{E,g) 



1. 



(2.102) 



In leading order in g, the corresponding relation may be derived on the basis of instanton calculus [see equations 
(7.66) and (7.70)]. 



2.3 Summary of the Quantization Conditions 

We summarize the conjectured quantization conditions. The general convention for the potential is indicated in 
equation (2. 10a). The A and B-functions are determined by equation (3.57), once the potential is fixed. All quan- 
tization conditions discussed here have been derived, in leading order in g, by instanton calculus. 

{Double-well potential.) For the double- well potential (2.10b), 



we have (2.21), 



H ■ 



2 

Bd„{E,g) 



exp 



AdUE,g) 



■■ £1 . 



(2.103) 



(2.104) 



The functions A^^ and B^iw defined in (2.23), and the perturbative quantization condition, from which the 
perturbative expansion of an eigenvalue may be obtained, is given in (2.25). 

(General symmetric potential.) We now consider a general symmetric potential with degenerate minima of the 
form (2.54), 

V{q) = ^q^ + 0{q^), V {q) = ^ (q - qo)^ + O {{q - qo)') , (2.105) 
where we have used (2.70). The quantization condition is (2.61), 



1 



B.y„(B,s) 



exp 



^sym 



{E,9) 



■ £1 . 



By contrast, the purely perturbative quantization condition is (2.62). The constant C is given in (2.58), 

/ 1 1 1 ' 



C 



Qo-q 



(2.106) 



(2.107) 
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{Asymmetric wells, degenerate minima.) For general asymmetric wells with degenerate minima [see equa- 
tion (2.63)], we have 



V{q) = \q^ + 0{q^) , V{q) = ^ a;' (? - qof + 0{{q- qof) . 
The generalized Bohr-Sommerfeld quantization formula takes the form (2.64), 

1 



(2.108) 



ra-B^{E,g))ra-B2iE,g)) 



1 

+ 2^ 



2a 

9 



Bi(E,g) 



20; 



B2{E,g) 



The perturbative quantization condition is (2.66). Here, we have [see equation (5.3)] 



exp ■ 



2uj 

1+U) 



go 



dq 



^/2V{q) Q w(go - q) 



(2.109) 



(2.110) 



{Periodic cosine potential.) For the periodic cosine potential V{q) = ^ (1 — cos 4(7), the conjecture then takes the 
form 



-Bpc(E,g) 



Apo(B,s)/2 



e p 
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Bpc{E,g) 



e p 



Apo(B,s)/2 



2cos</5 

T[\ + B^,{E,g)] ~ ^ ■ 



(2.111) 



The functions A^^ and B^^ are defined in (2.81), and the perturbative quantization condition can be found in (2.82). 

{Self-adjoint extension of the anharmonic oscillator with 0{v)-symmetry, at negative coupling.) For the 0{u)- 
anharmonic oscillator with the radial Hamiltonian (2.84), 



Hi = Hi{g) = -- 



1 + - r + or . 

2 r dr 2 ^ 2 



(2.112) 



2 Vdr, 

For g <0, the radial Hamiltonian being endowed with a self-adjoint extension, the quantization condition is (2.92) 

1 



r + J - B^iE,g,j))] r - .7 - B,{E,g,j))] 



+ 



BAE.gJ) 



exp{A,(E,g)) 
2tt 



= 0. 



(2.113) 



The functions and B^ are defined in (2.89). 

{Resonances of the anharmonic oscillator with 0{u)- symmetry, at negative coupling.) For g' < 0, the resonances 
of the radial Hamiltonian are conjectured to foUow the quantization condition (2.102), 



le 



-AAE,g) 



B^{E,g) 



exp 



. j + l + B,{E,g) 
1 TT :z 



r[| {j + l~B,{E,g))] 



1, 



T[^{j + l + B,{E,g))] 

{Fokker-Planck Hamiltonian.) An important special case of (2.92) is the Fokker-Planck Hamiltonian (2.94) 



(2.114) 



H 



FP 
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2d^2 ■ 2g 
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The corresponding conjecture for the quantization condition is given by the equation (2.96), 

1 f 2\^-^(^'«) exp(-AFp(£;,5)) 



r {-\B^p{E,g)) r (1 - \BF^{E,g)) 
The functions A^v and B-pv are defined in (2.99). 
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= 0. 



(2.115) 



(2.116) 
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Chapter 3 



Perturbative and WKB Expansions from 
Schrodinger Equations 



3.1 Orientation 

In the simplest examples these conjectures, based on the semi-classical evaluation of path integrals and instanton 
calculus, have obtained independent heuristic confirmations by considerations based on the Schrodinger equations, 
as we now explain. 

We write the Schrodinger equation as 

Hi,{q) = -f + -V{q)i,{q) = E^{q). (3.1) 

^ 9 

We consider only potentials V that are entire functions. This allows extending the Schrodinger equation and thus 
its solutions to the q complex plane. 

Moreover, we assume that V{q) has an absolute minimum, for convenience located at g = 0, where V{q) = 
h<l^+0{q^). 



3.2 General Considerations 

3.2.1 The Transition to Quantum Mechanics 

We recall in this chapter how the transition from classical to quantum mechanics can be understood in terms of 
matter waves whose wave fronts are given approximately by the 2-dimensional manifolds of constant classical 
action. The classical action Aq for a particle whose Lagrangian is not explicitly time-dependent, reads 

Aoit,r) = [ dt'L{r{t'),r{t')). (3.2) 
Jo 

One immediately verifies the relations 

dAo dAo dAo 

^='^' ^ = ^+^^' ^^°=^' ^'-'^ 
where p is the momentum of the classical particle. This can be rewritten as 

^ = -{pr-L) = -H, (3.4) 
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where H is the Hamiltonian. Using H =p'^ / (2m) + V, where V is the potential and m is the mass of the particle, 
the Hamilton- Jacobi equation is recovered, 

The momentum p = V^o is thus perpendicular to the manifolds of constant classical action Aq = constant. That 
is to say, the classically allowed trajectories are perpendicular to the surfaces of constant action. On the other hand, 
the wave fronts of a traveling wave are characterized by a constant phase, and the wave moves perpendicular to the 
surfaces of constant phase. Therefore, choosing the ansatz 

4>o{t, r) = exp(i^o(i, r)) (3.6) 

for the wave function, we are led to the conclusion that Aq should be constant on the surfaces of constant ^. 
Because ^ has to have dimension one, but Aq has the dimension of an action, the identification 

^o{t,r) = — (3.7) 

is suggested. If the dynamics of the matter wave were purely governed by the laws of classical mechanics, gener- 
alized to matter waves, then the equation (3.5) would adequately describe the time evolution of the wave packet 
described by the matter wave function tpQ, and the transition from classical to quantum mechanics would be com- 
pletely equivalent to the transition from ray to wave optics. In order to verify how well (3.5) is applicable to the 
description of a quantal wave packet, we choose the ansatz 

ilj{t,r) = exp(^^A{t,r)^ (3.8) 

for the quantum mechanical wave function and insert this into the Schrodinger equation. The result is 

This is not identical to (3.5). The third term on the right-hand side characterizes a diffusive process in the context 
of the Brownian motion of the quantal particle [22] and is responsible for the spreading of quantum mechanical 
wave packets. It may be interpreted as a quantum mechanical correction to the classical equations of motion of 
matter waves. The correction is of order h. An expansion in powers of ?l is a semiclassical expansion, or WKB 
expansion [23]. 

For a stationary state, the ansatz 

A{t,r)=a{r)-Et (3.10) 

leads to 

E=^^ + Vir)-^V^a. (3.11) 

We now specialize this equation to the one-dimensional Hamiltonian given by equation (2. 10) where, upon multi- 
plication by a overall factor g, the coupling g lakes the formal role of h with m = 1. We therefore have to replace 
E ^ g E, as well ash^g, and we also set a{q) = i J S{q') dq'. Then, 

gE=-ls\q)+V{q) + ^S'{q). (3.12) 

3.2.2 Riccati Equation 

According to the previous chapter, a convenient way to generate semiclassical expansions is to use the Riccati 
equation derived from the Schrodinger equation (3.1) by setting 

%) = -5W- (3.13) 
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The function S then satisfies (3.12) 

gS'{q)-S\q) + 2V{q)-2gE = 0. (3.14) 

Following the discussion of the appendix A, we also introduce two independent solutions '01,2 of equation (3.1) 
such that 

Mq) -t 0' Mq) 0, (3.15) 

and the decomposition (A. 13) with 

1(1 -I)- 

[The derivations of a few simple identities needed in this chapter are recalled in appendix A, but with different 
normalizations.] 

In the leading WKB expansion, the two wave functions ipi and ip2, for large q, i.e. in regions where V{q) > 
gE, may be expressed as exp(— 1 /g J dq y/2 V{q) — 2gE) and exp(l/ g J dq ^y2V{q) — 2gE), respectively (the 
positive square root being implied). 

Then, equation (3.14) translates into [see (A. 15)] 

gS'_ - Si - Si +2Viq) -2gE = 0, (3.17a) 
gS'+-2S+S- =0. (3.17b) 

We note that in equations (3.17) a change {g,E) i-^ {—g,—E) can be formally compensated by the change 
S-) I— > —S-). That is to say, S+ is even under {g, E) i— > {—g, —E), whereas S- is odd. 

It follows that in the sense of a series expansion 

S±{q,-g,-E) = ±S±{q,g,E). (3.18) 

The relations (3.17) allow to write the wave function in terms of S+ only, 

- I ri 



i,{q) = iS+)-'/' 



exp 



dq'S+{q') 



(3.19) 



3.2.3 Spectral Equation and Fredholm Determinant 

We now consider the resolvent operator 

R{E) = [H - E]-'^ . (3.20) 

The poles in E of its trace 

G{E) =TrR{E) = Tt[H - E]-\ (3.21) 

give the complete spectrum of the Hamiltonian H. Note that this expression may require some regularization to 
deal with the divergence of the sum over eigenvalues. The resulting ambiguity corresponds to adding a polynomial 
in E and does not affect poles and residues, i.e. it leaves the spectrum invariant. 

The trace of the resolvent G{E) is equal to the negative logarithmic derivative of the Fredholm determinant 

V{E) = det{H - E), 

= (3.22) 

in terms of which the spectral equation reads 

V{E) = det{H -E)^0. (3.23) 



This can be seen as follows. 



^lnl^(i^) = ^lndet{H-E) 



-Trj^ = -G{E). (3.24) 
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3.2.4 Spectral Equation and Logarithmic Derivative of the Wave Function 

In the normalization 3.1, the Schrodinger equation can be written as, 

(-d^ + 2 V{q)lg^ - 2E/g) ^{q) = , (3.25) 

which, in the notation of appendix A, implies L = 2 (H— E) / q, u(q) = 2 V{q) /g'^, z = 2 E/g and thus a different 
normalization for the resolvent (3.20). The corresponding diagonal elements r{q) of the resolvent [H—E]~^ satisfy 

lg2rr"-^ r'2 = 2 {V{q) -Eg)r^-1. (3.26) 

The equation defines r{q) up to a sign that is fixed by the behaviour for E — oo. The trace of R{E) is then given 
by 

G{E) = Tr —1- = [dq r{q), (3.27) 



H-E 

an expression that requires some regularization to deal with the large q divergences. 

Returning to the formalism of appendix A, one finds r{q) oc l/S+{q). A comparison between the different nor- 
malizations leads to 

r{q) = l/S+{q). (3.28) 

We then conclude that 

GiE) = J dq^^y \nV{E) ^-J dqS+{q). (3.29) 
The integral may diverge for |g| ^ oo and has then to be regularized. 

One way of writing the spectral equation, or quantization condition, is then [see equation (A.24)] 

2i^ }^+ j ^' ~ ~ ^' ^^''^ = + I ' 



(3.30) 



for > 0. In the case of analytic potentials, the contour of integration can be locally deformed in the q complex 
plane, and the limit e — > can then be taken. The spectral equation in terms of the logarithmic derivative of the wave 
function (3.13) becomes 

B{E,g) = -^ (f dzS+{z,g,E) = N+^, (3.31) 
2mg Jc 

where N has also the interpretation of the number of real zeros of the eigenfunction, and C is a contour that 
encloses them. 



3.2.5 Perturbative Expansion 

The purpose of the current chapter is to investigate the condition (3.31) in detail. Before embarking on this en- 
deavour, we remark that the rather elegant formulation discussed here is restricted, however, to one dimension 
and analytic potentials, which represents integrable systems from the point of view of classical mechanics. Our 
treatment bypasses the difficulties generally associated with turning points. This is the form that is useful here in 
the context of perturbative and WKB expansions. Recently, it has been pointed out by Rosenfelder [24] that non- 
analytic energy shifts of the form exp(— a/g) also occur in potentials which cannot be represented as polynomials. 

In general, canonical considerations which lead to the Bohr-Sommerfeld quantization conditions (see, e.g., chapter 
III of [23]), are replaced here by a simpler rigorous approach, based on the Riccati equation, that uses only the 
uniqueness of the wave function at g = 0. We thus assume that the potential has a unique minimum at g = 0. 
In order to see the connection with the uniformity of the wave function, we note that the condition (3.31) may be 
rewritten as 
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At (7 — > 0, this implies that the coefficient of g"^ in the Laurent expansion of S+ {q) is of the form —gN/ q where 
N is an integer. Idem est, we have In il){q) N Inq for g — > 0, a relation which is uniform under a sign change of 
q only if N is an integer. 

The expansion for g ^ alE fixed of the equations (3.17), or (3.26), leads to the expansion in powers of g of the 
function B{E, g) that enters the left-hand side of (3.31). This function enters in expressions like (2.22) and other 
quantization conditions. We first indicate how the function B{E, g) can be calculated (/) as a series expansion in 
powers of g bX E fixed, which leads to the perturbative expansion, and ( ii) at Eg fixed, which leads to the WKB 
expansion. We also discuss the relation between the two expansions. 

In view of (3.17), (3.18) and (3.19), we conclude that 

B{E,g) = -B{-E-g). (3.33) 

The latter reflection symmetry is purely formal — it is valid in the sense of power series — since it corresponds to 
change the quantum number N m—1 — N, which is unphysical. 

The perturbative expansion is obtained by expanding the solution of Riccati's equation (3.14) in powers of g (at E 
fixed): 

S{q) = Y.9''sk{q). (3.34) 

Setting 

U{q) = ^/2Vy), U{q) = q + 0{q'') , (3.35) 

we first find 

sM^Uiq), Siiq)^^^{U'{q)-2E). (3.36) 

The square root in (3.35) is to be understood such that the resulting function is uniform at the two minima of 
V, and positive in the region between the minima. For example, in the case of the double-well potential V{q) = 
q^ (1 - qf/2, we have U{q) = q{l- q). 

Higher orders are obtained from the recursion relation 

Sk{q) = ^ [s'^-M - > : ■sk-i{q)8i{q) | . (3.37) 

For example, 

U" 

S2(9) = ^-^77T + 7T77s -TTTTs- (3.38) 

Thus 

^+ = ^('^^ ' ^W) +5'(:^775-i^-:^) + (3-39) 

The contour integral (3.31) reduces to the residue at q = (where the potential has its unique absolute minimum). 
Thus, it depends only on the expansion of U{q) al q — Q. To find the perturbative spectrum, it is sufficient to 
expand the recursion relation (3.37) in powers of q. For two (or more) degenerate minima, this procedure gives rise 
to perturbative functions Bi,B2, ■ ■ ■ in the sense of chapter 2 [see also equation (3.57) below]. 

For instance, at order g, one finds B{E, g) = E. The next orders, for a general polynomial potential, are given in 
equation (B.U). The condition (3.31) then imphes that the eigenfunctions are uniform at g = 0, as they should be. 

There is an interesting connection to Borel summability. Note that even in situations where the perturbative ex- 
pansion of eigenvalues is Borel summable, it is not clear whether the function B{E, g) is Borel summable in g at 
E fixed. Indeed, the wave function ^(g) is unambiguously defined only when the quantization condition is satis- 
fied with N a non-negative integer. When E is not an eigenvalue, the solution of the Schrodinger equation is an 
undefined hnear combination of two particular solutions. 



-1(9)- 
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3.2.6 WKB Expansion 

As explained in chapter 3.2.1, the WKB expansion [25] is an expansion for 5 ^ at Eg fixed, in contrast to the 
perturbative expansion where E is fixed. In [26, chapter 6], it has been stressed that the loop expansion in quantum 
field theory is actually an expansion in powers of h ^ g, i.e. a generalized WKB expansion. We thus expand 
equation (3.14) in the complex g-plane in powers of g, at Eg fixed, starting from 

S{q) = So{q) , S^{q) = 2 V{q) -2gE. (3.40) 

For E > 0, the function 5*0 has two branch points g_ < g+ on the real axis. We put the cut between q- and q+ 
and choose the determination of 5*0 to be positive for q > q+. This ensures the decrease of wave functions on the 
real axis for \q\ — > 00, at least in the WKB approximation. 

Then [see equation (B.13)], 

S{q) = Y,9''Sk{q). (3.41a) 



fe>0 



One first obtains 



For fc > 1, one finds the recursion relation 



5,^1-^ (3.41b) 



fe-1 



Skiq) = isUiq) - j . (3.41c) 



At order g'^, for example, 



In the semi-classical limit, the contour C in the expression (3.31) encloses the cut of /S'o(g) which joins the two 
solutions of So{q) = (classical turning points). 

Note that equation (3.31) can also be written as 



exp 



1 

9 Jc 



. (3.43) 



The contribution of 5*1 as well as that of all Sj with odd j to the symmetric component 5*+ vanishes. This is 
evident as we inspect the symmetry condition (3.18), the properties of the leading approximation (3.40) and take 
into account the additional odd-power-of-gf prefactor that prevails in all Sj with odd j. Therefore, we may safely 
ignore all odd-order WKB approximants in the sequel; this simplifies the analysis carried out in appendix F. 

We now consider the WKB expansion of the resolvent. In the WKB limit Eg = 0{1), the solution of equation 
(3.26) at leading order, in the notation (3.14), is 



r{q) - 77^ = ^ G{E) = [ (3.44) 

So{q) ^2V{q)-2Eg J ^2V{q)-2Eg 



Formally, 



\nV{E) = ^J dqSoiq) = i y" dq ^2V{q)-2Eg, (3.45) 

but the divergence of the integral for |g| — > 00 has to be eliminated by subtracting a ^^-independent infinite 

constant. 

At order 5^, one finds [see equation (B.16)] 

S+{q) = So{q) + ff' - jir) + ^(9% (3.46) 
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and thus, after an integration by parts, 

\nV{E) = -J dg5o(g) + j dqS'^''{q) S^\q) + 0{g''). (3.47) 

In the present formulation, it is easy to relate the WKB and the perturbative expansions. The terms in the WKB 
expansion of the quantization condition (3.31) are contour integrals that, clearly, are regular functions of E at 
£■ = 0. To obtain the perturbative expansion from the WKB expansion, one thus expands the functions 5*2^ {q) in 
powers of E and calculates the residues. 

For example, at leading WKB order [see also appendix B.2 and equation (B.18)], 

One has then to evaluate the residue at the origin of 

-T^ i dqU^-^"" . (3.49) 
2i7r Jc 

More generally, replacing 5+ by its WKB expansion and expanding each term in a power series of Eg, one obtains 
the perturbative expansion of the function B{E, g). The WKB expansion corresponds to successive summations, 
to all orders in g, of the terms of highest degree in E of the perturbative expansion. A few terms are calculated in 
appendix B. 



3.3 Resolvent and Degenerate Minima in Specific Cases 
3.3.1 Potentials with Degenerate Minima 

We now consider potentials of double-well type, with two (not necessarily symmetric) degenerate minima. The 

symmetric case is displayed in figure 3.1, whereas the discussion here is actually applicable also to the case of 
minima with unequal curvatures [see equation (2.63)]. For E small enough, the function So{q) has four branch 
points g'l , . . . , (74 on the real axis. 

One expects, on intuitive grounds, to find two perturbative functions Bi{E,g) and B2{E,g) obtained by integrat- 
ing S+ around the two cuts [qi, (J2] and [q^, (J4] of figure 3.1. 

However, if one thinks in global terms, an obvious difficulty arises: a function So{q) with the corresponding two 
cuts is no longer an appropriate solution because it increases both for g ± 00, and thus the wave function cannot 
decrease both for q ± oc . Moreover, the contour integral of So caimot give the correct answer. For example, if the 

potential is symmetric, the integral vanishes. 

For symmetric potentials, the correct perturbative answer is obtained from the difference between the two integrals 
when the additional zero at the symmetry point for odd eigenfunctions is taken into account. One then finds 

2B{E,g) = l + 2[N/2], (3.50) 

where [N/2] is the integer part of N/2, and this is the correct perturbative answer. 

Moreover, quite generally the WKB expansion of the spectral condition has a singularity when 2E reaches the 
relative maximum of the potential between the degenerate minima, where two complex singularities of the So{q) 
pinch the real axis and the number of real singularities passes from two to four. 

Therefore, let us examine what happens if one starts from E large enough, where So has only two branch points 
on the real axis and no special difficulty arises, and proceeds by analytic continuation of the spectral condition. At 
leading order, before continuation, the function B{E, g) can be written as 

B{E,g)=B{E,g)=— dq ^2 Eg - 2V{q) + 0{g) , (3.51) 

'^9 Jq-(E) 
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where q± are the two real zeros of Eg — V. To avoid the singularity when E reaches the local maximum of the 
potential, we take g, and thus Eg slightly complex. Eventually, the initial contour is deformed into a sum C of 
three contours around the cuts: 



B{E,g) = Bi{E,g) + B2{E,g) + Bi2{E,g) 



with 



Bi2iE,g) = — / dq^2Ei 
ng ' 



«3 



2V{q)+Oig), 



(3.52) 



(3.53) 



where the cuts now are [qi , q2] and [ga , (74] . The function B12 is imaginary and the determination of the square root 
in the integral depends on the continuation. The integral can also be written as 



Bi2iE,g) = ^(f dq ^2 Eg - 2V {q) + 0{g) , 



(3.54) 



where C surround the function cut along [q2, q^]. Note that Bi2{E, g) as a function of E has a branch point at 
£J = 0. In a rotation about an angle 2-k around E = Q, B12 becomes B12 ± (Bi + B2). 
The quantization condition (3.31) or alternatively (3.43), 



N - 



1 

2 ' 



exp 



-- i dzS+{z) 
9 Jc 



1 = 0, 



(3.55) 



then becomes 

exp[2i7r {Br{E,g)+B2{E,g) + B^2{E,g))] + l = Q. (3.56) 

The contribution B12 contains the effect of barrier penetration and we choose the sign in front of it, which depends 
on the analytic continuation, to obtain decreasing corrections to perturbation theory. Finally, the contribution B\ + 
B2 can be absorbed into in a redefinition of -612 by going to the proper sheet in the E Riemann surface. 
Comparing with the instanton result, one infers the decomposition [in the notation of equation (2.64)] 

1 



9 Jc 



dzS+{z) = A{E,g)+ln{2'n:) 



.X:{lnr(l-B,(^,,))+B.(^,,)ln(-^)} 



(3.57) 



where the constants Cj are adjusted in such a way that A{E, g) has no term of order g^ and the coefficient of 5° 
is proportional to the expansion of B\ + B2. This equation, adapted to specific cases, almost directly leads to the 
conjectures discussed in chapter 2 and sunomarized in chapter 2.3 when combined with the quantization condition 
(3.55). 

The expansion for Eg small of its WKB expansion yields the function A{E, g). To identify with the WKB expan- 
sion, the function r(i — B) has to be replaced by its asymptotic expansion for B large: 

-lnr(i-B) ~ B\ti{-B)-B ~ B\n{-E) + --- (3.58) 

IT g ^0 



Therefore, at leading order. 



1 r ^ 

-Re i dzS+{z)= A{E,g) + \n{2'K)+y2Bi{E,g) 
3 Jc- i=i 



In 



yE_ 
2Ci 



1 



(3.59) 



Here, we have assumed that the determinations of all functions of E are such that A{E, g) is real in perturbation 

theory, something that is possible since the ambiguities are proportional to 2\'k{Bi + B2): in the WKB expansion, 
the discontinuity comes from the asymptotic expansion (3.58) of the sum of the F functions, which also has a dis- 
continuity — 2i7r(-Bi +B2), showing the consistency of the whole scheme. Moreover, the contribution proportional 
to ln(— g) leads to the combination hi{gBi), which is globally even in E, g as it should be [because it results from 
an integral over 5"+ with the symmetry property (3.18)]. 

A calculation of A{E, g) at a finite order in g requires the WKB expansion and the asymptotic expansion of the 
F-function only to a finite order. For example, the expansions up to order (2.24b) and (2.81b) of the function 
A{E, g) for the double-well and cosine potentials, which have initially been determined in part by numerical 
calculations, are reproduced by the expansion of the two first WKB orders (the calculations for the double-well 
potential are given in appendix B). 
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3.3.2 General Symmetric Potentials 

We first consider general symmetric potentials: V{—q) = V{q). As is evident from the discussion in chapter 2.1 
above and chapter 4.2.1 below, it is necessary to introduce a further quantum number in this particular case, 
which is the parity. States are classified according to the principal quantum number N and the parity e. We again 
combine the formalism of chapter 3.2.6 and appendix A. Using equation (A.28), and taking care of the different 
normalization, one infers 



Ga{E) = TvP[H-E]-' =2 / ^-^exp — / dy S, 

Jo '^+[1) L 9 Jo 



(y) 



(3.60) 



where P is the reflection operator that acts here according to Pipiq) = V'(~9) [cf- equation (4.11) where the 
reflection symmetry is about q = 1/2 in contrast to the case of a center of symmetry at g = discussed here]. 

An alternative integrated form is [equation (A.30)] 

V.iE) = exp i-Tr PHH - E)] « ^^,(0) = . (3.61) 

At leading order in the WKB limit, one finds —{2E)~^, the harmonic oscillator result forE — > — oo [equation 
(A.39)]. 

In the case of the symmetric double-well potentials, this situation again is more subtle, because one has to consider 
two different WKB approximations near the two wells, number in this case A strategy consists in calculating [27] 

Ga{\{E+,N + i^-,iv)) ~ ^ . (3.62) 

We thus consider the Schrodinger equation in the normalization (3.1) with E = N -\-^ and call qo> Q the position 

of one of the minima of the potential: V{q) ^ \{q — qoY- 

For g > 0, at leading order in the semi-classical limit, 

S+{q) = ^m{q) - 2Eg = U{q) - + 0{g^) , (3.63) 
where U{q) is defined as [see equation (3.35)] 



C/(g) = V2V%), U{q)=q + 0{q^). (3.64) 
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We choose the branch cut of the square-root function such that U{q) becomes a function regular at g = go and 
positive for g > go- The approximation (3.63) is valid for |g — go| large. 

The solution V'jv(g) can be written as 



V'iv(g) 



We introduce the definitions 



Qo {qo - q) 
X exp 



q-qo 



U{q) 



cxp 



1 

.9 Jo 
1 

) q' -Qo 



dg' U{q') 



rgo 

4 / dqU{q), 
Jo 



C = 



ql exp 



90 



dg 



U{q) q + qo q - qo 



(3.65) 

(3.66) 
(3.67) 



When comparing to (5.6), it is necessary to remember that we are dealing with a symmetric potential V{—q) = 
V{q) in the current section, and we integrate from the center of symmetry at g = to one of the minima of the 
potential. By contrast, in (2.60) and in (5.6), the integration extends from one of the minima to the other. Also, the 
convention for the square root is the opposite U = —U for g inside the region between the two minima. 

Based on (A.27a), we can write ra(g) as 



ra(g) = C 



_ „-N-l/2 a/2g 2 



with 



(3.68) 
(3.69) 



??7v(g)~(go-«)^e-(«-«°) 

an expression valid for 1 <C |g — go| <C y^. 

The leading contribution to the integral comes from the neighbourhood of g = go, where the harmonic oscillator 
approximation is relevant. The corresponding eigenfunction with norm ^/n is 

1 . ...r _„^/, 2^/2 



VN{q) 



q—yoo J TV! 



2^/2 vM 

Therefore, r?jv(g) has to be approximated by 

7?iv(a) - ViV!(ff/2)^/VAr((g - go)/v^). 

We conclude 

G^{\{E^^^E^J)) ~ 2C-''-'/^m{2/gfe''/^'^ Jdq<p%{{q-qo)/V9) 
~ 20F3 (7-^-1/2 TV! (5/2)^ e«/29. 

Using equation (3.62), one obtains the one-instanton contribution to the energy difference 



(3.70) 



(3.71) 



2ttNI 



2C 



N+l/2 



-a/2g 



(3.72) 



(3.73) 



a result consistent with the expansion of the zeros of expression (5.18). 

Note that to go beyond leading order by this method is simple for N fixed, but not for generic N because the WKB 
expansion yields the perturbative eigenfunction under the form of an expansion only vahd for large arguments. 

To compare with instanton calculations, we calculate the contribution of the twisted resolvent at leading order. 



TrPln(ff-£;) = X! ( 



En - SEn/2 - E 
En + 6En/2- E 



-E 



SE 



■JV 



JV 



En — E 



(3.74) 
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where E has been assumed to be not too close to an eigenvalue. At leading we can replace Epf hy N + |. Thus, 

TrPln{H-E) 

O 1 /or'X -^+1/2 poo 



1/2 



—j e-^/^s j d/3 exp [P{E - i) + 2 Ce-^^/g] . (3.75) 
We note that if we evaluate the integral in the hmit g — > 0_, we obtain 

TrPln(i/-i;) ~ -2— ^ exp [i7r(^- i)] ( — j ~ E), (3.76) 

an expression that appears in equation (4.60) and has the correct poles and residues. Of course, the real part, 

e-a/29 /2C\^ 

TrPln(iJ-j;) ~ -2 cos(7r(^ ^ 1/2)) — r(i - £;), (3.77) 

v27r \ 9 J 

has the same poles with the same residues, but has zeros for E <Q. 

3.3.3 (9(z/)-Symmetric Hamiltonians in the Radial Coordinate 

A general O (i/)-symmetric Hamiltonian, after diagonalization of the angular momentum, decomposes into a set of 
Hamiltonians Hi, where I is the angular momentum, of the form 



/ay i^- 1 d l{l + v-2) 
\ dr / r dr 



-V{r). (3.78) 
.9 



We assume that V{r) is an analytic function of r^, in such a way that the Schrodinger equation can be continued to 
the complete real axis. Eigenfunctions are then odd or even, and they behave hke r' for r — > 0. The eigenfunction 
ipi^N corresponding to angular momentum I and quantum number N thus has I + 2A'' zeros on the real axis. The 
quantization condition (3.31) is replaced by 

— i dz-^^^=Z + 27V, (3.79) 

2m Jc Wi,n{z) 

where the contour C now encloses all zeros of the eigenfunction tpir) on the real axis (this includes r < 0). 
The transformation 

^ = -V)-^, (3.80) 
ipyr) g 2r 

leads to the Riccati equation 

gS'{r) - S\r) + 2V{r) - 2Eg + g^ "/^^ = (3.81) 

withj = i//2 + l-l. 

Introducing the two functions S±, we can write the quantization condition as 



/ dzS+iz)=j + 2N + l. 

^17^5 Jc 



(3.82) 



In the perturbative expansion of the 0{v)-symmetric Hamiltonian, the effect of angular momentum appears only 
at order g. In the parameterization (B.6) [see also equations (7.48), (B.ll) and (B.12)], this implies the following 
relation between Bv{E, g, j) and B^wiE, g): 

B4E,g,j)=2Bi^(^^,-g^+^a2 (f-l) g + 0{g^). (3.83) 
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The WKB expansion of the 0{u)-symmetric Hamiltonian is now an expansion sXEgmA {j^ — 1/4) fixed. The 
leading order is 

S{r) = ^U^{r)-2Eg+^ 0' " i) ' ^^'^"^^ 
Again, the perturbation series is recovered by simply expanding in powers of g. 
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Chapter 4 

Instantons in the Double- Well Problem 



4.1 Orientation 

The conjectures presented in chapter 2 have been initially motivated by semi-classical evaluations of path integrals, 
more precisely by sunnmation of leading order multi-instanton contributions. 

Therefore, we now explain, first in the example of the double- well potential, how multi-instanton contributions can 
be summed explicitly, yielding results consistent with the conjectured form (2.22) of the secular equations. The 
reason for still emphasizing instanton calculus in situations now rather well understood by other methods, is that 
the arguments can be generaUzed to other examples where our understanding is more Umited. 

4.2 General Considerations 
4.2.1 Partition Function and Resolvent 

The path integral formalism allows calculating directly the quantum partition function, which for Hamiltonians 
with discrete spectrum has the expansion 



N=0 

Eigenvalues are labeled En (N e No), with Eq being the ground-state energy and En < En+i - The trace G{E) 
of the resolvent of H is related to the partition function by Laplace transformation: 



where we have to assume initially that E < Eq (with Eq being the ground state energy), to ensure the convergence 
of the integral for /3 — > oo. 

The poles of G{E) then yield the spectrum of the Hamiltonian H. From G{E) one can also derive the Fredholm 
determinant X'(-B) = det{H — E), which vanishes on the spectrum (chapter 3.2.3). 

Note that the integral (4.2) does not necessarily converge at /3 = because the eigenvalues may not increase fast 
enough, as the example of the harmonic oscillator shows: 



which diverges. However, G{E) — G(0) is defined, and fixing the irrelevant constant term we can set [see equation 




(4.1) 




(4.2) 




(4.3) 



(A.32)] 



Go,,.{E) = -xl,{\-E) 



d_ 
dE 



Inr(i-iJ) . 



(4.4) 
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Note that for the double-well potential, one can separate eigenvalues corresponding to symmetric and antisymmet- 
ric eigenfunctions by considering the two functions 

oo 

Z±{P) = Tr [i (1 ± P) e-^^] = ^ e-'^^i-" , (4.5) 

JV=0 

where P is the parity operator (2.11). The eigenvalues are then poles of the Laplace transforms (e = ±): 

G,{E) = f d/3e^^Z,(/3). (4.6) 
Jo 

4.2.2 Path Integrals and Spectra of Hamiltonians 



In the path integral formulation of quantum mechanics, the partition function is given by 

[dq(t)\exp 

lq{-[i/2)=q(fi/2) 



Z{p)<x [ [dq{t)]exp 

Jn(-6/2)=n(B/2) 



■S{q{t)) 
9 



(4.7) 



where the symbol /[d(;(t)] means summation over all paths which satisfy the boundary conditions, that is closed 
paths, and S{q) is the Euclidean action: 



/•/3/2 

<S(g)= / [\e{t) + V{q{t))]cit. (4.8) 



In [28], it has been stressed that the formulation in terms of the EucUdean action follows naturally from an analytic 
continuation of the usual definition of the path integral to imaginary time. 

In the case of the symmetric potential (2.10) which has degenerate minima, it is actually necessary to also consider 
the quantity 

^a(/3) = Tr {Pe-"") cx / [dq{t)] exp ^ 



--S{q(t)) 



(4.9) 



/g(-/3/2)+g(/9/2) = l 

where P is the parity operator (2.11). Then, eigenvalues corresponding to symmetric and antisymmetric eigen- 
functions can be derived from the combinations 

^±(/3)=Tr [i(l±P)e-^^] =i(Z(/3)±Za(/3)). (4.10) 

This construction is not necessary (even impossible) in the case of asymmetric wells discussed below in chapter 5. 
Because there is no parity in the asymmetric case, the quantity Z^{j3) does not find a natural representation in 
terms of eigenfunctions in contrast to the equations (4.13) and (4.11) below, and it is only the quantity Z{(5) as 
defined in (4.7) which contributes to the partition function. Of course, Z{(}) is composed exclusively of paths in 
which the starting and the endpoints are identical, and one-instanton effects therefore do not contribute to Z{P). 

We adopt here a convention, in which the proportionality sign in the notation (4.7) impUes that the final integration 

over Q = q{(3/2) is understood. 

Because the only irreducible representations of the parity operator are one-dimensional, the nondegenerate eigen- 
functions have definite parity, 

P(j),,N{l~q)^e^eM{q). (4.11) 

Here, £ = ±1 is the eigenvalue of the parity operator defined according to (2.1 1). It is instructive to consider the 
eigenfunction decomposition of the spectral partition function Z{p), 

Z{(3) = fdq Y,<Pe,N{q)<P:,N{Q) exp(-/3£;,,jv) • (4.12) 

Here, we assume a Hamiltonian whose spectrum is discrete (no continuous spectrum). The quantity Za{P) then 
finds a representation as 

2a{P) = Yl I dg0e,iv(g)0:,iv(l-9) eM-(iEe,N). (4.13) 

e,N ^ 
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Now, Zj^ (/3) is determined by the even eigenfunctions alone, 

oo „ 

= ^ / dg.^+,jv(g)<^;,;v('Z)exp(-/3^+,jv) 

(4.14) 



Af=0' 

2 

whereas Z- (/3) is determined by the odd eigenfunctions. 



1-^ n 

Z_{P) = ^ / dg0_,jv(g)r-,iv(5)exp(-/3i;-,iv) 
Z{f3)-Z,{(3) 



(4.15) 



4.2.3 Perturbation Theory 

Perturbative expansions, that is expansions in powers of 5 for 5 0, can be obtained by applying the steepest 
descent method to the path integral. Saddle points are solutions qc{t) to the Euclidean equations of motion (Eu- 
clidean equations differ from the normal equations of classical mechanics by the sign in front of the potential). 
When the potential has a unique minimum, located for example at g = 0, the leading saddle point is qc{t) = 0. A 
systematic expansion around the saddle point then leads to a purely perturbative expansion of the eigenvalues of 
the Hamiltonian of the form 

00 

E{g) = Y,Eig^. (4.16) 

1=0 

In the case of a potential with degenerate minima, one must sum over several saddle points: to each saddle point 
corresponds an eigenvalue and thus several eigenvalues are degenerate at leading order. Because the potential 
(2. 10) is symmetric, the lowest eigenvalue, and more generally all eigenvalues which remain finite when g goes to 
zero, are twice degenerate to all orders in perturbation theory [see equation (2.5)]: 

00 

E±M9) « E^n\9) ^ E . (4.17) 



4.2.4 Instanton Configurations 

Eigenvalues which remain finite when g goes to zero can be extracted from the large /3 expansion. In the infinite /3 
limit, leading contributions to the path integral come from paths which are solutions of the Euclidean equations of 
motion that have a finite action. In the case of the path integral representation of Zi^ifl), constant solutions of the 
equation of motion do not satisfy the boundary conditions. Finite action solutions necessarily correspond to paths 
which cormect the two minima of the potential (see figure 4.1). 

In the example of the double-well potential (2.10), such solutions are 

q,{t) = (1 + e±(*-*°)) , S{q,) = 1/6 . (4.18) 

Since the two solutions depend on an integration constant to, one finds two one-parameter families of degenerate 
saddle points. 

Non-constant (in time) solutions with finite action are called instanton solutions. Since the main contribution to 
the action comes from the region around t = to, one calls to the position of the instanton. 

The corresponding contribution to the path integral is proportional, at leading order in g and for (3 ^ 00, to 

g-i/ (6g) jjjyj, jj, non-perturbative. It is also proportional to /? because one has to sum over all degenerate saddle 
points and the integration constant to varies in [0, /?] for /3 large but finite. The complete calculation involves taking 
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the time to as a collective coordinate, and integrating over the remaining fluctuations in the Gaussian limit. One 
finds that the two lowest eigenvalues are given by {s = ±) 



l_e-i/69(i + 0(^)). 



(4.19) 



4.3 Multi-Instantons 



Taking into account E'^^^ (g) and E'^^^ (g), one obtains for the functions (4.5) an expansion of the form 

ZM ^^e-/^^'-^^^') ~ e-^^' Y: ^ {^o T- (4.20) 

Thus, the existence of a one-instanton contribution to eigenvalues impUes the existence of n-instanton contributions 
to the functions (4.5), proportional to /3". 

For /3 finite, the path integrals indeed have other saddle points which correspond to oscillations in the well of the 
potential —V{q). In the infinite (3 limit, the solutions with n oscillations have an action n x 1/6 and thus give 
contributions to the path integral of the expected form. 

However, there is a subtlety: naively one would expect these configurations to give a contribution of order /? (for /? 
large) because a given classical trajectory depends only on one time integration constant. This has to be contrasted 
with the expansion (4.20) where the n^^ term is of order 

Indeed, one discovers that the Gaussian integration near the saddle point involves the determinant of an operator 
that has eigenvalues which vanish exponentially in the large- /3 limit. The divergence has the following origin: in the 
large-/? Umit, the classical solution decomposes into a succession of largely separated instantons and fluctuations 
which tend to change the distances between instantons induce an infinitesimal variation of the action. It follows 
that, to properly study the hmit, one has to introduce additional collective coordinates that parameterize all config- 
urations close to solutions of the Euclidean equation of motion, even though they have a sightly different action. It 
is then easy to understand where in the expansion (4.20) the factor /3" originates from: Although a given classical 
trajectory can only generate a factor /?, these new configurations depend on n independent collective coordinates 
over which one has to integrate. 

To summarize: we know that n-instanton contributions do exist. However, these contributions do not correspond, 
in general, to solutions of the classical equation of motion. They correspond to configurations of largely sepa- 
rated instantons connected in a way which we shall discuss, which become solutions of the equation of motion 
only asymptotically, in the limit of infinite separation. These configurations depend on n times more coUective 
coordinates than the one-instanton configuration. 
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4.4 Specific Calculations 
4.4.1 Instanton Interaction 

We now briefly explain, still with the example of the double-well potential, how multi-instanton contributions to 
the path integral can be evaluated at leading order for g 0. 

In the infinite /3 Umit, the instanton solutions can be written as 

q±{t) = /(T(t-io)), (4.21a) 
fit) = l/(l + e*) =l-/(-i), (4.21b) 

where the integration constant to characterizes the instanton position. 

We first construct the two-instanton configuration. Actually it is convenient to now call instanton a solution which 
goes from to 1 and anti-instanton a solution which goes from 1 to 0. Then, the relevant configurations are 
instanton-anti-instanton pairs. These configurations depend on one additional time parameter, the separation be- 
tween instantons, decompose in the limit of infinite separation into two instantons and for large separation must 
minimize the variation of the action [29, 30]. For this purpose, we could introduce a constraint in the path inte- 
gral fixing the separation between instantons (see chapter C.3), and solve the equation of motion with a Lagrange 
multipUer for the constraint. Instead, we use a method which, at least at leading order, is simpler and more intuitive. 

^ ^ ^ 



-10 -5 5 10 15 

Figure 4.2: The two-instanton configuration. 



We consider a configuration [t) that is the sum of instantons separated by a distance 6, up to an additive constant 
adjusted in such a way as to satisfy the boundary conditions (figure 4.2): 



qS) = fit - e/2) + fi-t - e/2) -i = fit- e/2) - fit + e/2) , 



(4.22) 



where fit) is the function (4.21) (and of course all configurations deduced by time translation). This path has the 
following properties: it is continuous and differentiable and when 6 is large it differs, near each instanton, from 
the instanton solution only by exponentially small terms of order e~^. Although the calculation of the correspond- 
ing action is straightforward, we give here some details to show that the ansatz (4.22) applies to more general 
symmetric potentials. 

It is convenient to introduce some additional notation: 



uit) = fit -e/2), 
vit) = uit + e) , 

and thus qc = u — v. The action corresponding to the path (4.22) can be written as 



(4.23) 



= J dt[^ql + Viq,)] 
= 2x^ + j dt [-iiv + Vi; 



U — V) 



Viu) - Viv)] . 



(4.24) 
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The parity of Qc allows us to restrict the integration to the region t>0, where v is at least of order e After an 
integration by parts of the term vii, one finds 



1 



= - + 2 |t>(0) m(0) + y dt [vu + V{u-v) -V{u) -V{v)]'i . (4.25) 

One then expands the integrand in powers of v. Since the leading correction to S is of order e~^, one needs the 
expansion only up to order v^. The term Unear in v vanishes as a consequence of the u-equation of motion. One 
obtains ^ 

S{qc) - ^ ~ 2 v{0) w(0) + 2 1^ °°dt[^v^ V"{u) - | ^"(0) v^] | . (4.26) 

The function v decreases exponentially away from the origin so the main contributions to the integral come from 
the neighbourhood oft = 0, where u=l + 0{e-^/'^ and thus V"{u) ~ V"{1) = V"{0). Therefore, at leading 
order the two terms in the integral cancel. At leading order, 

v{0) u(0) e"^ (4.27) 

and thus 

S{q,) = ^-2e-' + 0{e-^') . (4.28) 

It will become clearer later why the classical action is needed only up to order e~^. In analogy with the partition 
function of a classical gas (instantons being identified with particles), one calls the quantity — 2e~^ the interaction 
potential between instantons. 

Actually, it is simple to extend the result to P large but finite. Symmetry between 9 and P — 9 then implies 

S{qc) = \- 2e"^ - 2e"('^"^) + negUgible contributions. (4.29) 
This expression can be verified by calculating its extremum as a function of 6. One obtains 

9, = 13/2, ^ S{q,) = i - 4e-'^/2 + O (e"'^) . (4.30) 

For the same Hamiltonian, one can calculate for (3 finite the action corresponding to a solution of the equation of 
motion with one oscillation. One finds 

5(gc) = ^-2e-^ + 0(e-2/3) . (4.31) 

Both results are consistent. Indeed, to compare them one has to replace (3 by [3/2 in equation (4.31) and multiply 
the action by a factor 2, since the action corresponds to a trajectory described twice in the total time (3. 

We now examine the variation of the action due to the multi-instanton configuration. Specifically, we show that if 
we modify infinitesimally (for 9 large) the configuration to further decrease the variation of the action, the change 
r{t) of the path is of order and the variation of the action of order e~^^ at least. Setting 

q{t) = q,{t) + T{t) (4.32) 

and expanding the action up to second order in r{t), one finds 

S{q, + r) = S{q,) + J [q,{t) r{t) + V' {q,{t)) r{t)] dt 

+ 1 J dt [f^t) + V"{q,) T^t)] + O {[r{t)f) . (4.33) 

In the term linear in r(t), one integrates by parts r{t), in order to use the property that qc{t) approximately satisfies 
the equation of motion. In the term quadratic in r(t), one replaces V" by 1, since r{t) is expected to be large only 
far from the instantons. One then verifies that the term Unear in r is of order while the quadratic term is of 
order 1. A shift of r to ehminate the linear term would then give a negligible contribution, of order e~^^. 



42 



CHAPTER 4: INSTANTONS IN THE DOUBLE-WELL PROBLEM 



4.4 Specific Calculations 



We now consider an n-instanton configuration, i.e. a succession of n instantons (more precisely, alternatively 
instantons and anti-instantons) separated by times 9i with 

Y.1=A = P- (4.34) 

At leading order, we need only consider "interactions" between nearest neighbour instantons. Other interactions are 
negligible because they are of higher order in c^^. This is an essential simplifying feature of quantum mechanics 
compared to quantum field theory. The classical action Sc{9i) can then be directly inferred from expression (4.29): 

n 

S^iOi) = I -2j2e-'' + O (e-(^*+^^)) . (4.35) 

i—l 

Note that for n even, the n-instanton configurations contribute to Tre"'^^, while for n odd they contribute to 
Tr (Pe~^^) (P is the parity operator). But all contributes to the combination (4.5). 

(Remark.) Since we keep in the action all terms of order e~^, we expect to find the contributions not only to the 
two lowest energies but also to all energies which remain finite when g goes to zero. 



4.4.2 The n-Instanton Contribution 



We have calculated the n-instanton action. We now evaluate, at leading order, the contribution to the path integral 
of the neighbourhood of the n-instanton configuration [7, 31]. We expand the action up to second order in the 
deviation from the classical path. Although the path is not a solution of the equation of motion, it has been chosen 
it in such a way that the linear terms in the expansion can be neglected. The Gaussian integration involves then the 
determinant of the second derivative of the action at the classical path 



M{t', t) = 



+ V"{qc{t)) 



S{t-t'). 



(4.36) 



The operator M has the form of a Hamiltonian with a potential that consists of n wells asymptotically identical to 
the well arising in the one-instanton problem, and which are largely separated. At leading order the corresponding 
spectrum is, therefore, the spectrum arising in the one-instanton problem n-times degenerate. Corrections are ex- 
ponentially small in the separation. Simultaneously, by introducing n collective time variables, we have suppressed 
n times the zero eigenvalue and generated the Jacobian of the one-instanton case to the power n. Therefore, the 
n-instanton contribution to the combination (4.5) 



Z,(/3) = iTV [(l + £P)e-'5^] 



(e = ±1), can be written as 



l/6s 



n V 



Jei>o ^ ' , 



i=l 



(4.37) 



(4.38) 



All factors have already been explained, except the factor /3, which comes from the integration over a global time 
translation, and the factor 1 /n, which arises because the configuration is invariant under a cyclic permutation of the 
di. Finally, the normalization factor cT^^I"^ corresponds to the partition function of the harmonic oscillator. Odd-n 
instanton effects contribute positively to Z^"^ and negatively to zi"' (/3). 

If the instanton interactions are neglected (but instanton configuration with an arbitrarily large number of tunnel- 
ings are allowed), the integration over the ^j's is straightforward and the sum of the leading order n-instanton 
contributions 



g) « g) = e'^^/^ + ^ Zj") g) 



can be calculated: 



Ze(/?,<?)«^Se(/?,5) = e- 



-/3/2 



n 



n=l 



-1/63 



/ (n- 1)! 



(4.39) 



(4.40) 
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with 

EM = \ + 0{g)- ^e-V65 (i + o (g)) . (4.4I) 

We recognize the perturbative and one-instanton contribution, at leading order, to -Eg^o (5) > the ground state and the 
first excited state energies. 

In view of problematic issues related to the instanton interaction, the summation of these effects and necessary 
analytic continuations which have to be performed in order to calculate eigenvalues, we add here a small discus- 
sion. To go beyond the one-instanton approximation, it is necessary to take into account the interaction between 
instantons (see chapters 4.4.1 and 4.4.2). Unfortunately, if one examines expression (4.38), one discovers that 
the interaction between instantons is attractive. Therefore, for g small, the dominant contributions to the integral 
come from configurations in which the instantons are close. For such configurations, the concept of instanton is no 
longer meaningful, since the configurations cannot be distinguished from fluctuations around the constant or the 
one-instanton solution. 

We should have expected such a difficulty. Indeed, the large-order behaviour analysis has shown that the perturba- 
tive expansion in the case of potentials with degenerate minima is not Borel summable. An ambiguity is expected 
at the two-instanton order But if the perturbative expansion is ambiguous at the two-instanton order, contributions 
of the same order or even smaller are ill-defined. To proceed any further, we must first give a meaning to the sum 
of the perturbative expansion. 

In the example of the double-well potential, it is possible to show that the perturbation series is Borel summable for 
g negative, by relating it to the perturbative expansion of the 0(2) anharmonic oscillator. Therefore, we define the 
sum of the perturbation series as the analytic continuation of this Borel sum from g negative to 51 = 1^1 ± iO. This 
corresponds in the Borel transformation to an integration above or below the real positive axis. We then note that, 
for g negative, the interaction between instantons is repulsive and, simultaneously, the expression (4.38) becomes 
meaningful. Therefore, we first calculate, for g small and negative, both the sum of the perturbation series and 
the instanton contributions, and perform an analytic continuation to g positive of all quantities consistently. In the 
same way, the perturbative expansion around each multi-instanton configuration is also non-Borel summable and 
the sum is defined by the same procedure. 

At the same time, it should be remembered that the various analytic continuations and cancellations of imaginary 
parts, while illustrating the internal consistency of the approach, are not really necessary in order to obtain the 
energy eigenvalue via generalized resummation of the resurgent expansion: it would suffice to say that the energy 
eigenvalue is obtained by assigning to every divergent series in g, i.e. to every perturbative expansion about the 
n-instanton contributions, the real part of the Borel sum associated with that (factorially divergent) series. The 
real part of the Borel sum is also obtained by evaluating the Laplace-Borel integral using the principal-value pre- 
scription. Although intuitively clear, we would like to remark that this principal- value prescription — of course — 
corresponds to the integration contour Cq as outiined in [17]. 

While, on the one hand, the necessity for an analytic continuation may seem a little obscure at first sight (see for 
example remarks on p. 241 of [32]), this procedure is, on the other hand, a mathematically well defined concept. 
In the context of the double-well problem it gives a defined meaning to each of the factorially divergent, nonalter- 
nating power series occurring in equation (2.13). In the end, given the fact that even mathematical proof [10, 1 1] 
may be questioned by inquisitive minds, the soundness of the concept must be judged by its predictive power. Ex- 
empli gratia, one may ask the question whether the energy levels predicted by the analytically continued instanton 
expansion are able to accurately reproduce numericaUy determined levels. This is indeed the case, as shown in 
chapter 8 below. 

4.4.3 The Sum of the Leading-Order Instanton Contributions 

Let us now introduce two useful parameters [cf. equations (2.27) and (2.28)] 

/x=--=exp(x(5)) , A=^e-V65, (4.42) 
9 si lix 
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in such a way that the "fugacity" ^{g) of the instanton gas, which is half the one-instanton contribution at leading 
order, can be written as [see (2.27)] 



(4.43) 



where the quantity ^(9) also occurs in a form multiplied by the parity e, and a sign change in the determination of 
-y/— /i is therefore equivalent to a change e — e [see also equation (4.55) below]. 

As discussed before, we now assume that initially g is negative. The Laplace transform 



f-OO 

Jo 

of the n-instanton contribution (4.38): 

/3e-/3/2 (-ieA^M)' 



/ 5(^0, d^.exp 



.5 7=1 



(4.44) 



(4.45) 



yields the leading contribution Gi"'\E) to the trace Ge{E) of the resolvent [equation (4.6)]. In order to factorize 
the integral over the 9i, we introduce a complex contour integral representation for the 5-function, 



27ri 



i=l 



In terms of the function 



I{s,iJ,) 



p+oc 

Jo 



exp (s9- ide~^) dO , 



(4.46) 



(4.47) 



•Zi"^ (/?) can be rewritten as 



/Je-/5/2 (-ieA)" 



27rin 



100 



(4.48) 



In view of (4.44), we have 



POO 

G^J'Xe) = / d/3e^-^zi")(/3) 
Jo 



d/3 



^e/3(£-l/2) (_i£A)" 



Jo 27rm 

ioo 

dE J 27rin 

—ioo 

ioo 

J 27rin 



100 



JO 



8 + 1/2- E 



dE n 



(4.49) 



Alternatively, one may observe that in (4.44), the integral over (3 is immediate and the integrals over the 9i then 
factorize. One obtains 



G^I^\E)r.ti^^lX{E-\,^r . 



(4.50) 
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The function T{s, fi) is the only function needed for these leading-order calculations in one dimension. A first way 
to calculate it, is to expand 

p + oo °° ( i\N ( T\N ,,N+l/2 

^M=l E^/^""^^^e^^-"'^d.= E^^. (4.51) 
"'u jv=o ■ JV=0 

The expansion defines a real meromorphic function of s with poles at non-negative integers. The poles and residues 

can be trusted because they depend only on the large 9 behaviour. However, after analytic continuation to g > 0, 
the function grows as e"'' = e^/^, a behaviour due to the integration near = that cannot be trusted and cannot 
be correct. 

Let us now evaluate the integral (4.47) in the Umit /x — > -|-oo, and thus 51 — > 0_. We change variables, setting 
/X =t, and the integral becomes 

pH r+oo 

I{s,i^) = dii-i-^e"* = /i*+^/M dtt-^-'e-' +0 {c-^'/y/JI) . (4.52) 

Jo Jo 

We thus obtain {s = E -1/2) 

I{s, m) « /x^+i/^ n-s) =(^-^y r(^^-Ey (4.53) 

a meromorphic function with the same poles and residues as the initial expression. For /i +00, the difference is 
exponentially small and both expressions are equivalent, but after analytic continuation to g > 0, the asymptotic 
form (4.53) has now an acceptable behaviour. Therefore, our ansatz is that the estimate (4.53) gives the correct 
leading behaviour of the true function. 

In view of (4.50), the generating function GeiE, g) of the leading-order multi-instanton contributions (4.39) then 
is given by [we use = — -|- x)] 

g,{E,g) = Y,Gt\E) = -—\nA,{E) (4.54) 

with ^ 

A,{E) = l + siXl(^E-^.f?j « 1 + (^-^^ rQ-i?^ , (4.55) 

where we now explicitly see that a change in the determination of ^— /x is equivalent to a change £ 1— > — e. Because 
QeiE, g) approximates the trace of the resolvent. 



g,{E,g)^T:r[^^—^] , (4.56) 



we have 



Qe{E,g) ^ -^TrMH,-E)] 
= -^\n[det{H,-E)] 

« -^lnA,(^). (4.57) 

Here, -ff^ denotes the Hamiltonian restricted on the Hilbert space spanned by the eigenvectors of H with parity 
e. Therefore, the functions A.s{E) w det(ile — E) are directiy the sums of the leading order multi-instanton 
contributions to the products ^^{E) of even or odd eigenvalues, 

^e{E)<x]l(l-^^ , (4.58) 
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but, unlike T>s{E), they are meromorphic functions in E because the zero-instanton contribution has not yet been 
included at all; an expansion of (4.55) in powers of A starts with a term of order 0(X). Correcting for this effect, 
we add to Ge{E, g) the trace of the resolvent of the harmonic oscillator, and we thus divide ^e{E) by r(i — E). 
This amounts to the replacement 

ge{E,g) ^ g,{E,g) + Go,c.{E) 

- -4 in (4.59) 

dE r{^-E) ^ ^ 

This cancels the poles, and (E) becomes an entire function: 

The first term is simply the Fredholm determinant (A. 33) corresponding to the harmonic oscillator. We note here 

that the sum of all instanton contributions, in the leading approximation, i.e. neglecting subleading terms in the in- 
stanton interaction, simply yields some sort of one-instanton (but i?-dependent) correction to the spectral equation. 

Since A is small, zeros of the equation As{E) = are close to eigenvalues of the harmonic oscillator: 

Ee,N = N+^+0{X), N>0. (4.61) 
The zeros of the function (4.60) can then be expanded in a power series in A: 

E,M9)=J2^N\9){-sXr. (4.62) 

n 

One obtains from a unique equation the multi-instanton contributions to all energy eigenvalues Eg^j^{g) of the 
double-well potential at leading order in g, that is to say all the coefficients e N,nk[) in the notation of the equations 
(2.13) and (2.15). Concrete results for the double-well potential are presented in chapter 8.7. 

The appearance of a factor In (/ in equation (2.39) can now be simply understood by noting that the interaction 
terms are only relevant for g"^ of order 1, that is 6 of order — \ng. 

Assuming that an equation of the form (4.60) holds beyond leading order, we conclude that the argument of the 
F-function must be such that the equation reproduces at zero instanton order the perturbative expansion: E has to 
be replaced by the function B{E, g). If we further assume that for B large there is some connection with the WKB 
expansion, we infer that {—2/g)^ must also be replaced by {—2/g)^ to reconstruct a term of the form B ln{Eg). 
This is the origin of the conjecture (2.21). 

This concludes the first part of the treatise of multi-instanton effects in quantum mechanics. In the second part [33], 
we intend to generaUze the instanton calculations to a wider family of potentials. 
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A.l Inverse Schrodinger Operator: Matrix Elements 

In this appendix, we recall here the derivation of a few classical results that we have used throughout this article. 
We first consider the hermitian positive differential operator 

L = -dl+u{x), (A.1) 

where u{x) has the form of a potential, which remains strictly positive for \x\ — > oo. 

We recall the derivation of the differential equation satisfied by the diagonal matrix elements of its inverse R = 

Matrix elements R{x, y) of R satisfy the Schrodinger equation 

{-dl + u{x)) R{x, y) = 8{x - y) . (A.2) 
We recall that R{x, y) can be expressed in terms of two independent solutions of the homogeneous equation 

{-dl + u{x))ip{x)=Q. (A.3) 
We denote hy ip\,'^2 two solutions, partially normaUzed by 

and, moreover, satisfying the boundary conditions 

(^i(a;) — > for a; — > — oo, y2(a;) — > for a; — > +oo . (A.5) 
Then, it is easily verified that R{x, y) is given by 

R{x,y) = fi{y) f2{x) 0{x - y) + ipY{x) Lp2(y) 0{y - x) , (A.6) 
where 9{x) is the usual Heaviside step function. The diagonal matrix elements 

r{x) = R{x,x) = ipi{x) ip2ix) (A.7) 

satisfy 

r'{x) = ip[ {x) (P2 (x) + (pi{x) (P2{x) , 

r"{x) = 2{^[{x)^'2{x)+u{x)r{x)) , (A.8a) 
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where we have used the equation (A.3). 

One then verifies that, as the consequence of equation (A.4), r{x) satisfies the non-linear differential equation 

2r{x) r" (x) - r'^ (x) - 4 u{x) (x) + 1 = . 



A quantity of special interest is the trace G of the R: 

^ ~ J da;r(a;). 

A.2 Riccati's Equation 



We now set 



in equation (A.3). We find 



(7{X) 



ip{x) ' 



cr'(x) -cr^(x) +u{x) = 0. 

To the two solutions 1^1,2 correspond two functions (Ti 2- We introduce 

Cr+ = 5 (fT2 - (Tl), fT_ = i((T2 + CTi). 

The Wronskian condition (A.4) implies 

1 1 



a+(x) 



2 (fii{x) (fi2{x) 2r{x) 



The equation satisfied by ai^2 implies 



ct'_ - al - a'i + u{x) = , 

C7+ - 2 C7+ cr_ = . 



The second allows expressing (t_ in terms of a+ and, thus, 

1 



V2 



Vi = 



exp 



exp 



(A.9) 



(A.10) 



(A. 11) 
(A. 12) 

(A. 13) 
(A. 14) 



(A.15a) 
(A. 15b) 



(A. 16a) 
(A.16b) 



A.3 Resolvent and Spectrum 

We now substitute u (a;) u{x) — z, setL = iJ — 2; and apply these results. The operator 
the resolvent of the Hamiltonian H. Riccati's equation becomes 



a'ix) - C7^(a;) + u{x) -z = 0. 



Then, 



da' 2 ^'^ 1 

dz dz 



One first useful consequence is the identity 

9(7+ (a; 



dz +M^)M^)=g^ 



1 da-ix) 



2cr+(x) dz 



L ^ (z) becomes 
(A. 17) 

(A.18) 
(A.19) 
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Then, 



9(7+ (a;) 



G{z) = 'TrL-\z) = j Axr{x) =- j dx ■ (A.20) 



Integrating, one infers 

/ \ \ +00 



In'D(^) = lndeti(2;) = j da;cr+(a;) = 2 



(A.21) 



The function V^z) is the Fredholm determinant of the operator H ~ z: Assuming a discrete spectrum, we have 

V{z) oc Jl (1 - z/zn). (A.22) 

N=0 

Note that the infinite product may require some convergence factor, and the expressions in equation (A.21), in 
general, require some large-a; regularization. 

The spectrum of H is given by the solutions of the equation ^^{z) = 0. An alternative and useful way of writing 
this spectral equation is 

lim liiV(zN - is) - InDf^jv + ie) = iTT(2N + 1) , (A.23) 

and, therefore, 

-— lim / dx[(7{x, zn — ie) — (T+{x, zn + is)] = N + ^. (A. 24) 

2?,7r £^0+ J 

Finally , the large x behaviour then determines the solutions completely and one finds 

^"^^•"^ - ' ^^^dy^liy), (A.25a) 

dy^liy) (A.25b) 



dz ifiKx) 

1 r 



dz ifli 



{General symmetric potentials.) We now assume u{x) = u{—x). We then choose 

(pii-x) = (p2ix) = ip{x). (A.26) 
Instanton calculus suggests that it is then also interesting to consider the quantities 

ra{x) = R{x,-x)=ip'^{x)e{x)+(p^{-x)e{-x), (A.27a) 
Ga{z) = Tr PL-^{z) = J dxra{x) = 2 dxip^ix), (A.27b) 

where P is the a; i-> —x reflection operator. 

Using equations (A. 16a) and (A.4) taken at a; = to normalize, we obtain a first expression 

f°° dx \ r 

Ga{z)= — TT^^P ~^ dya+{y) . 

Jo cr+W L ^0 

We then combine with equation (A.25) to obtain 



(A.28) 



In particular, integrating over z, one finds an expression for the generalization of the Fredholm determinant: 

exp [-TrPln(if - .)] = TT f ^^^"l oc = ^ = -^4^, (A.30) 



an expression that can be easily understood since for even eigenfunctions <^'(0) vanishes, while ip{0) vanishes for 
odd eigenfunctions. 
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A.4 Example: Harmonic Oscillator 

As an illustration, let us apply this formalism to the harmonic oscillator with u{x) = and z = 2E. 
The trace of the resolvent for the harmonic oscillator formally reads 

^osc.(^)=E ^+V2-i^ ' (A.31) 

JV=0 ' 

which diverges. However, the first derivative of G{E) is defined. Integrating one can choose 

GosciE) = -i^il - E), (A.32) 

where il>{z) is the logarithmic derivative of the F-function. Finally, from Gosc. {E) one derives the Fredholm deter- 
minant 

V,,,XE) = l/r{^-E), (A.33) 
which vanishes on the poles of the F-function. 

The solution of the Schrodinger equation that decreases for x +00 is proportional to 

where the contour encloses the real negative axis. The asymptotic form for a; ^ —00 is given by steepest descent: 

ip{x,E) ~ -^{-2x)-''-^/^ smn{E+he''"/\ (A.35) 

X— >-oo -^TT 

For X +00, the integral is dominated by the neighbourhood of the origin: 

V{X,E) ~ a:^-y^—}-^e-^'/\ (A.36) 



Then, 



Moreover, 



V(E) (X rr^ oc . (A.37) 

^ ' r(E+i) m-E) 



<p(0, E) = ^ 2-^-1/2 sin[7r(i? + i)] F Q - ^E^ , (A.38a) 

^'(0, E) = ~ sin[7r(£ + i)] F (^^ - ^E^ , (A.38b) 

and thus the expected result is obtained. 

For example, in the case of the harmonic oscillator, 

GaiE) = ^ - ^ ^ ^) = + i + Oil/E% (A.39) 
an expansion that is reproduced by the WKB approximation. 
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B.l Perturbative Expansion 

For illustration purpose, we first calculate a few terms of the perturbative expansion of the function B{E, g), using 
the expansion (3.37). We then verify the relation between perturbative and WKB expansions. 

We expand the Riccati equation (3. 14) in powers of g. 



k>0 



Setting 
we first find 



U{q) = ^/27(g) , 
1 



so{q) = U{q), si{q) = 



2U{q) 



{U'{q)-2E). 



Higher orders are given by solving the recursion relation 

1 



fc-i 



At order g"^, one finds 



sk{q) = 



S2{q) = 



2U{q) 



s'k-iil) -^Sk-iiq)si{q) . 



{2U"U - 3U'^ + 4EU' - 4£;2) 



8f/3(g) 

We then simply need the residues at g = 0. 

As we show later, it is often convenient to parameterize the potential as the solution of the equation 



(B.l) 

(B.2) 
(B.3) 

(B.4) 
(B.5) 

(B.6) 



Beyond perturbation theory, this implies simple well or symmetric double well potentials. For example, the quartic 
double-well potential corresponds to /? = 1 — 4u and the cosine potential to p = 1 — 4u^. 

Up to the order of g^, the solution has the expansion 



V{q) 



1 

2' 
+ 



1 



■ ai q 



1 



1 



ai Q!2 



■ as 5 



q'. 



(B.7) 
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This suggests an alternative method of calculation. Introducing the parameterization 

S = a+iU)+gU'a-{U)=U + 0{g), 
one obtains another form of Riccati's equation: 

(T^(7i) - 2cr+(u)cr_(u) = 0, 

[p{u){a'_{u) - a^_{u))+^p'{u) ct_(w)] - ct^(w) + u'^-2gE = 0. 



(One may first calculate cr^.) Then 



B{E,g)=- 



2i-Kg 



du 



(7+[u) 



and the residue at g = becomes the residue at u = 0. 

In terms of the parameters (B.6), one obtains up to the order g'^, 

1 3 



B{E,g) ^ E 
1 

■- a4 -f- 



-a2 



16 



16 



64 



+ —T^a4, + — aia3 + —a 



15 



25 



16 



32 



64 



128 "1"^+ 1024"^ 



aUE 



(B.8) 



(B.9a) 
(B.9b) 



(B.IO) 



(B.ll) 



One immediately verifies that it agrees with the results obtained for various special potentials. In the case of an 
even potential (ai = as = 0) and finite angular momentum, the expression at the same order becomes [see also 
equation (3.83)] 



B{E,g) = E + a2 



1 



17 



1 



1 



3 



20 



40 



(B.12) 



B.2 WKB Expansion: Second-Order Calculation 

{General considerations.) The WKB expansion of B{E, g) is an expansion in powers of g at gE fixed. Therefore, 
we expand S{q) in powers of g at gE fixed (or alternatively equation (3.26)): 

S{q)=Y.9''Sk{.q). (B.13) 

Then, 

and the recursion relation becomes 

Sk{q) = (sl_,{q) -'j2Sk-l{q)Si{q)^ • (B.15) 

At order g^, one finds 

iS"' 3<S"^ 
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B.2 WKB Expansion: Second-Order Calculation 



For the order contribution to the functions A[E^ g) and B{E, g), we need only the integral of 52. Integrating 
by parts, one finds 

/'"'Si ^ 

{Perturbative and WKB expansions.) One recovers the contributions to the perturbative expansion by expanding 
Sk {q) in powers of Eg and calculating the residues at this minimum of the potential. 

We expand the first WKB term in powers of gE: 

So = Uy(2gEU-T /^"7r l^ • (B-18) 

t'o r(n + i)r(-i) 

One has thus to evaluate [see also the equations (F.9) and (F.38)]: 

-— f dqU'-^^ = -—d>du^=, (B.19) 
2i7r Jc 2i7r / ^/p{u) 

where the parameterization (B.6) has been introduced. 

For the next WKB order, we calculate the integrand of (B.17) more explicitly 

Expanding So, one finds [see also the equations (FlO) and (F.39)]: 



(B.20) 



One thus has to evaluate 



_ J_ / dgt/'2t/-2"-3 = --L /d«y^u-2"-3. (B.22) 
2in Jc 2 1 TT y 

(A special class of potentials.) We follow [6]. It is possible to calculate many terms of the WKB expansion for the 
class of potentials such that, in the parameterization (B.6), 

p{u) = 1 - 4u'" . (B.23) 

Here, m = 1 corresponds to the double-well potential, m = 2 to the cosine potential, for m = 3, 4 the potentials 
are meromorphic elliptic functions. Higher values no longer correspond to meromorphic functions (except for 
TO = 6 where is still meromorphic). The case m > 3 odd corresponds to potentials which diverge on the real 
axis, m even to periodic and finite potentials on the real axis. 

One verifies that the term of order g' in the expansion is of the form Pi{u) / u'^^~^ , where Pi is polynomial. The 
expansion of B thus involves only the integral 

L{iy) = --^ <f dqU^-^r (B.24) 

We change variables q i-^ u = U{q) in (B.24) and find 

L„(z/) = - J- I du{l- 4m™)-^/^ u^-^^ . (B.25) 
2i7r 

Expanding the root in powers of u, we obtain 

L„(,y) = -24(^-1)/™ ^ " ^ . for 2{u-l) = (mod m) . (B.26) 

r(i)r(i + 2ii^) 
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B.2 WKB Expansion: Second-Order Calculation 



For example, 

^'""Jc ^0 r(n + i)r(-i) 

~ ~ r(^ + 2)r(-i)r(fc + i)r(i) ' ^^'^^^ 

where means the sum over k, including only those terms for which A: to is even. 
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Appendix C 

Multi-Instantons 



C.l Determinant 



In this chapter, we give some indication about a few additional technical ingredients that are involved in the 
calculation of multi-instanton contributions. We can write the operator M defined by equation (4.36) as 



M 



in which v{t) is a potential localized around t = 0: 

v{t) = O (e-l*l) , \t\^oo 



and ti are the positions of the instantons. 
We want to calculate 



det M Mq-^ = det ■ 



Using the identity In det = Tr In, we expand the r.h.s. in powers of v{t): 



|l + [-(d/dt)Vl] 



In det MMq-^ = ^ 



(-1) 



fc=i 



k 



/n 



A(wi — U2) 



H = l 



with the definition 



A{t) = {0 -{d/dty + l 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



It is clear from the behaviour of v{t) and of A(t), that when the instantons are largely separated, only the terms in 
which one retains from each potential the same instanton contribution survive. Therefore, 



In det MMq"^ = 



n 



E 



(-1) 



k+l 



k=l 

for \ti — tA » 1 . 



Y\ dwj A(Mi — U2) v{u2) ■ ■ ■ A(Mfe — Ml) v{ui) 



(C.6) 



We recognize n times the logarithm of the one-instanton determinant. 



56 



APPENDIX C: MULTl-INSTANTONS 



C.2 Instanton Interaction 



C.2 Instanton Interaction 



We assume, as in chapter 5, that the potential has two degenerate minima at the points a; = and x = xq with 



V{x) = ^x^ + 0{x^) , 

V{x) = i _ ^^^)2 + ((^ _ ^^)3^ 

Let us write the one-instanton solution qc{t) which goes from to go = xq/ y/g as 

qc(t) = f{t)/^. 

We choose the function f{t) in such a way that it satisfies 



xo-fit) - 

fit) \/C^e* fort-^-oo. 

By solving the equation of motion, it is easy to calculate the constant C^j'. 



xluj^/^'+^^ 



exp 



2aj 
1+w 



f 

Jo 



dx 



1 

X 



uj{xo — x) 



(C.7a) 
(C.7b) 

(C.8) 

(C.9a) 
(C.9b) 

(C.IO) 



We recognize the constant (2.73). 

We now construct instanton-anti-instanton pair configurations q{t) which correspond to trajectories starting from, 

and returning to, (7 = go or = 0. Since we want also to consider the case of two successive instantons, we assume, 
but only in this last case, that V{x) is an even function and has therefore a third minimum at a; = — xq. 
According to the discussion of chapter 4, we can take as a two-instanton configuration 



q,{t) = ^{U{t)+£f-{t)), 



±1 



(C.ll) 



with 



f+{t)=f{t-ei2), f_(t)^f{-t-0/2). (C.12) 

Here, is a measure of the instanton separation. The case e = 1 corresponds to an instanton-anti-instanton pair 
starting from g = go at time —00, approaching g = at intermediate times and returning to go. The case e = —1 
corresponds to a sequence of two instantons going from — go to go. Finally, for the classical trajectory which goes, 
instead, from the origin to go and back, we can take 

Q2{t) = [fit + 9/2) + f{e/2 -t)- xo] /V5 • (C.13) 

We now calculate the classical action corresponding to gi (t). We separate the action into two parts, corresponding 
at leading order to the two instanton contributions: 



5(gi) =5+(gi)+5_(gi) 



(C.14) 



with 



Jo 



-,q! + -v{y^qi{t)) 

^ .9 



dt, 
dt. 



(C.15) 



The value t = of the separation point is somewhat arbitrary and can be replaced by any value which remains 
finite when Q becomes infinite. We then use the properties that for d large /+(t) is small for f < 0, and f-{t) is 
small for < > 0, to expand both terms. For example, for 5+ we find 



- / dt 
9 Jo 



lifim+viuit)) 



+e 
1 

+ 2 



Ut)U{t) + v'{uit))f.{t) 
f'_{t) + v"{u{t))f_{t) 



(C.16) 
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C.3 Multi-Instantons from Constraints 



Since /_ {t) decreases exponentially, only values of t small compared to 6/2 contribute to the last term of equation 
(C.16) which is proportional to V". For such values of t, we have 

(C.17) 

For the terms hnear in /_ (t), we integrate by parts the kinetic term and use the equation of motion 

m = v [fit)] . (C.18) 

Only the integrated term survives and yields 

/' + 00 



/ 

Jo 



dt 



f+it) + v'if+it)) f.{t) = -fi-e/2) fi-e/2) . (C.19) 



The contribution <S_ can be evaluated by exactly the same method. We note that the sum of the two contributions 
reconstructs twice the classical action a. We then find 



9 

with 



S{qi) = - 2a-2ef{-e/2)f{-6/2) + - 



(C.20) 



a= ^/2V{x)dx. (C.21) 
Jo 

Replacing, for 6 large, / by its asymptotic form (C.9), we finally obtain the classical action 

= [2a - 2Cee-^ + O {e'^^)] (C.22) 

and thus the instanton interaction. 

Following the same steps, we can calculate the classical action corresponding to q2 (t). The result is 

5(92) = ^ {2a - 2 [f{e/2) - xo] f{e/2) + •••}, (C.23) 
which for 6 large is equivalent to 

S{q2) = ^ [2a - 2{C/w)e-'^^] . (C.24) 

Finally, if we consider the case of a finite time interval f3 with periodic boundary conditions, we can combine both 
results to find the action of a periodic trajectory passing close to g = and q = qo- 

S{q) = g-^ [2a - 2C (6"''+^ + e-'^Vw)] , (C.25) 

in agreement with equations (5.5) and (5.8). 



C.3 Multi-Instantons from Constraints 



Although multi-instanton configurations do not correspond to solutions of the equation of motion, it is nevertheless 

possible to modify the classical action by introducing constraints and integrating over all possible constraints. The 
main problem with such a method is to find a system of constraints which are both theoretically reasonable, and 
convenient for practical calculations. 

One can, for instance, fix the positions of the instantons by introducing in the path integral (in the example of the 
double-well) 



1 



/n j^tql{t--ti) ^ j dtqs,{t-U){q{t)-q,,{t-ti)) 



dti 



(C.26) 



where ti are the instanton positions and a succession of ± indicating instantons and anti-instantons. One then 
uses an integral representation of the (5-functions, so that the path integral becomes 

2 \ " r 



\q+\ 



2m 



JJ dti d\i I [dq{t)\ JJ exp [-S{q, \i)] with 

1=1 i=l 



n ^ 

S{q, \i) = S{q) +Y.\i dt q,, {t - ti) {q{t) - g,, {t - ti)) . (C.27) 

i=l 
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Appendix D 



Dispersion Relations and (non-)Borel 
Summability 



D.l Dispersion Relations 

(The simplest example.) The nonaltemating {g > 0) factorially divergent series 

oo 



is generated by expanding in powers of g tlie following integral, 



m = [ 

Jo 



ooibi e 



di- 



1 



l-gt 



exp(— . 



(D.l) 



(D.2) 



The imaginary part due to the (half-)pole alt = l/gis 



Im/(ff) = ±- exp(-l/ff). (D.3) 
Because f{g) has a cut along the positive real axis, one may write a dispersion relation in an obvious way as 



= -27ri Res 



=20 Z — Zq 



+ 



f 

Jo 



dx 



X - Zo 



{f{x + ie)-f{x-ie)) 



{Dispersion relation.) The equality 



/(^o) = - / 



oo -j^ 

dx Im f(x) 

X- Zo 



(D.4) 



(D.5) 



follows immediately. Let us now start from this dispersion relation as the fundamental property of a function 
which has a cut along the positive real axis and fulfills sufficient conditions in order to ensure the convergence of 
the integral. Let us assume furthermore, that /(zo) can be expanded as a formal power series. 



f{zo) ~ ^ a„ zj . 

n=0 

Expanding the right-hand side of (D.5) into a power series in zq, 

' — Zf\ X \ X y 



X — Zo X V X . 

n=0 



(D.6) 



(D.7) 
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D.2 A Simple Example of non-Borel SummabiUty 



we immediately have 

1 r°° 1 

an = - dx—^lmf{x). (D.8) 

In many cases, possible convergence problems near a; = are eliminated due to a nonanalytic factor of the form 
exp(— 1/a;) that enters into Im f{x). 

(A calculational example.) We now replace x ^ g. For a function whose imaginary part is 



Im f{g) = i exp ( - — ) {l + bg + eg" + dg'' + ...), (D.9a) 
the evaluation of the integral on the right-hand side of (D.8) leads to 

an = n! 1 + - + + +...). (D.9b) 

TT \ n n[n + l) n[n + l)(n + 2) J 

This example shows that correction terms to the imaginary part of order g to the imaginary part along the positive 
real axis (the cut) correspond to subleading corrections to the factorial growth of the "perturbative coefficients" a„ 
of relative order n~^. The correction terms can be summarized in a natural way in terms of an inverse factorial 
series. This consideration is the basis for the "matching" discussed in chapter 8.2. 



Figure D. 1 : The four roots of equation (D. 14). 



D.2 A Simple Example of non-Borel Summability 

Let us try to illustrate the problem of non-Borel summability with the example of a simple integral, which shares 
some of the features of the problem in quantum mechanics which we have studied in chapter 4. We consider the 
function 

1, 



1 



+ 00 



(D.IO) 

9 

where V{x) is an entire function with an absolute minimum at a; = 0, V{Q) = 0. For g small I{g) can be calculated 
by steepest descent, expanding V around q = Q: 



l{g) = Y.lk9^. (D.ii) 



fc>0 

It is easy to write a finite dimensional integral of the form (D. 10) as a generalized Borel or Laplace transform. 



We integrate over q: 



I{g) = -^JdqdtS [V{q^) - t] e-*/^. (D.12) 



1 f°° 1 
Hg) = / dte-*/s y , , (D.13) 
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D.2 A Simple Example of non-Borel SummabiUty 



in which {xi{t)} are the solutions of the equation 

V[x,it)]=t. (D.14) 

When the function V{x) is monotonic both for x positive and negative, the equation (D.14) has two solutions for 
all values of t and the equation (D.13) is directly the Borel representation of the function I{g), which has a Borel 
summable power series expansion. 



I I 
I I 
i I 
I t 
I I 
I I 
i 



Figure D.2: The different contours in the a;-plane. 



We now assume, instead, that V{x) has a second local minimum which gives a negUgible contribution to I{g) for 
g small. A simple example is 



V{x) 



^x 



l.,s^r + a) + lx\ 



i < a < 1 , 



(D.15) 



which has a minimum at a; = 1. Between its two minima the potential V{x) has a maximum, located at x 
whose contribution dominates the large order behaviour of the expansion in powers of g: 



a. 



Ik oc 



r(fc) [Via)]-'' , V{a) > 



(D.16) 



(in the example (D.15) V{a) = (1 — o/2)/6) and the series is not Borel summable. 

The naive Borel transform of I{g) is obtained by retaining in equation (D.13) only the roots of equation (D.14) 
which exists for t small. The singularities of the Borel transform then correspond to the zeros of V"{x). 

For the potential (D.15) the expression (D.13) has the form 



1(9) 



1 



/2^ 



+00 



dte-*/9 



1 



e{V{a)-t)9{t-V{l)) 
\V'{xs{t))\ 



V'{x,{t))\ \V'{x,{t))\ 



\V'{x4t))\ 



(D.17) 



with the definitions (see figure D.l): xi{t) < < X2{t) < a < X3{t) < 1 < X4{t). 

The idea of the analytic continuation is to integrate each contribution up to t = +oo following a contour which 
passes below or above the cut along ihe positive real axis. This means that we consider X2 (t) to solve the equation: 



V[x2{t)]=t±ie. 



(D.l 8) 



The sign is arbitrary. Let us, for instance, choose the positive sign. We then have to subtract this additional contri- 
bution. We proceed in the same way for 0:3 (t) for t > V{a). Since X2 (t) and X3{t) meet att = V{a), the analytic 
continuation will correspond to take for xs{t) the other solution 



V[x3{t)]=tTi£- 



(D.19) 
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D.2 A Simple Example of non-Borel SummabiUty 



We thus have to subtract from the total expression the contributions of two roots of the equation. But it is easy to 
verify that this is just the contribution of the saddle point located at a; = a, which corresponds to a maximum of 
the potential. 

Therefore, we have succeeded in writing expression (D.17) as the sum of three saddle point contributions (see 
figure D.2). There is some arbitrariness in this decomposition which here corresponds to the choice e = ±1. 

In the complex x plane, we have replaced the initial contour C on the real positive axis, by a sum of three contours 
Ci, C2 and C3 corresponding to the three saddle points located at 0, a, 1. 
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Appendix E 



Degenerate Minima: Energy Splitting and 
Schrodinger Equation 



We consider the Hamiltonian 



H 



dl + V{x) 



.1^2 
2 



(E.1) 



where the potential is analytic on the real axis, even: V{x) = V{—x) and has an absolute minimum located at 
X = ±Xm (choosing Xm > 0) where it vanishes. Eigenfunctions and eigenvalues can be calculated in a pertur- 
bative expansion in each well of the potential. There are twice degenerate to all orders. We call tpN{x) the WKB 
eigenfunction corresponding to the expansion near the well around 



(E.2) 



Then the second eigenfunction with the same perturbative energy is iPn{—x) and it is a WKB approximation 
near —Xm- We know that the true eigenfunctions are even or odd. To calculate the energy difference SE^ = 
E_^_^N — E_^N we construct trial wave functions (j)N,e{x), which are even or odd depending on the value e = ±1, 



iPn{x) + etpNi-x), 



for a < X , 
for — a < X < a , 
for X < —a , 



(E.3) 



where Xm > a > 0. At the points x = ±a the function tpN{x) is exponentially small in the semi-classical limit 
andV'Ar(-Q:) < tpN{oe)- 

The trial wave functions are continuous, satisfy locally the Schrodinger equation, but are not differentiable at the 
points X ±a. We now calculate 



E 



{(t>N,e\(l>N,E) 



(E.4) 



Integrating by parts, one finds 

r+oo 



I 

J —I 



dx 



^{cl>'^Jx)Y + V{x)^%^,{x) 



= E 



N 



+ 00 



It follows, that up to smaller exponential corrections, 

V'Jv(a)V'jv(-") + V'Af(-a)V'jv('^) 



5E 



N 



Jdx-ip^ix) 



(E.5) 



(E.6) 
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where only the neighbourhood of x = Xm contributes to the integral in the denominator. We note that since 
iI)n{±x) are two eigenfunctions with the same energy, the quantity 



W{a) 



(E.7) 



is a Wronskian and thus independent of a. The calculation then follows closely the lines of chapter (3.3), where 
the inverse quantity is calculated. Equation (E.6) can easily be generaUzed to the Schrodinger equation in higher 
dimensions. 

We now apply this result to the Schrodinger equation in the normalization (2.10). At leading order in the semi- 
classical limit, the eigenfunction ^n{<i) can be written, for < q < q^, as 



■>pN{q) 



q"" Vq/UW) exp 



1 f 
- / dq'U{q') 
9 Jo 



X exp 



{N- 



dq' 



1 



1 



U{q') q' 

At this order ijj'N{q)/tpN{q) ~ —U{q)/g. The value of the Wronskian W{q) thus is 



(E.8) 



w{q) = -tpNiq)ipNiqo- q) 

= ^ (jN+l/2 ^-a/2g 

9 



i>N{qo-q) tpNiq) 



with [see equations (2.58), (2.60) and (3.35)] 

a = 2 dqU{q), 
Jo 

"/■««/ 1 
Jo 



C 



9o exp 



1 



(E.9) 

(E.lOa) 
(E.lOb) 



U{q) q qo-q, 

For the denominator, we need the perturbative eigenfunction near q = Q with the same normalization. We thus set 
q = x^ and obtain 

^n{x^9) -y/Ve-^'/'. (E.11) 

g^O 

This result, however, yields the coefficient of c^^^/^ in the WKB limit q = fixed, that is x large. Actu- 
ally, we know that at leading order the wave function is simply the eigenfunction of the harmonic oscillator. The 
eigenfunction with norm ^/^^ is 



iPn{x) = 



2N/2 



We conclude 
Therefore, 



Vjv(a;V5)'-ViV! (5/2)^/2 



dq^pUq) ^ Nl {g/2) 



N 



(E.12) 

(E.13) 
(E.14) 



Finally, we notice that in the q variable and with the form (E.l) of the Schrodinger operator the spectrum is E/g. 
The one-instanton contribution to the energy difference follows: 



5EN{g) 



2C\ 



JV+l/2 



-a/29 



(E.15) 



s/2^N\ \ g J 

a result consistent with the expansion of the zeros of expression (5.18). 

Note that to go beyond leading order by this method is simple for fixed, but not for generic N because the WKB 
expansion yields the perturbative eigenfunction under the form of an expansion only vahd for large arguments. 
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